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ABSTRACT 

Superwinds from starburst galaxies are multi-phase outflows that sweep up and 
incorporate ambient galactic disk and halo gas. The interaction of this denser material 
with the more diffuse hot wind gas is thought to give rise to the Ovi emission and ab- 
sorption in the far Ultraviolet (FUV) and the soft thermal X-ray emission observed in 
superwinds. In this paper we present high-resolution hydrodynamical models of warm 
ionized clouds embedded in a superwind, and compare the Ovi and soft X-ray prop- 
erties to the existing observational data. These models include thermal conduction, 
which we show plays an important role in shaping both the dynamics and radiative 
properties of the resulting wind/cloud interaction. Heat conduction stabilizes the cloud 
by inhibiting the growth of Kelvin-Hclmholtz and Rayleigh- Taylor instabilities, and 
also generates a shock wave at the cloud's surface that compresses the cloud. This 
dynamical behaviour influences the observable properties. We find that while OVI 
emission and absorption always arises in cloud material at the periphery of the cloud, 
most of the soft X-ray arises in the region between the wind bow shock and the cloud 
surface, and probes either wind or cloud material depending on the strength of con- 
duction and the relative abundances of the wind with respect to the cloud. In general 
only a small fraction (< 1%) of the wind mechanical energy intersecting a cloud is 
radiated away at UV and X-ray wavelengths, with more wind energy going into ac- 
celerating the cloud. Clouds in relatively slow cool winds radiate a larger fraction of 
their energy, which are inconsistent with observational constraints. Models with heat 
conduction at Spitzer-levels are found to produce observational properties closer to 
those observed in superwinds than models with no thermal conduction, in particular 
in terms of the Ovi-to- X-ray luminosity ratio, but cloud life times are uncomfortably 
short (< 1 Myr) compared to the dynamical ages of real winds. We experimented 
with reducing the thermal conductivity for one set of model parameters, and found 
that even when we reduced conduction by a factor of 25 that the simulations retained 
the beneficial hydrodynamical stability and low Ovi-to- X-ray luminosity ratio found 
in the Spitzer-level conductive models, while also having reduced evaporation rates. 
Although more work is required to simulate clouds for longer times and to investigate 
cloud acceleration and thermal conduction at sub-Spitzer levels in a wider range of 
models, we conclude that thermal conduction can no longer be ignored in superwinds. 

Key words: ISM: clouds - ISM: jets and outflows - galaxies: starburst - ultraviolet: 
galaxies - X-rays: galaxies 



1 INTRODUCTION 

Superwinds are multi-phase, loosely-collimated galaxy-sized 
outflows with measured velocities in excess of several hun- 
dred to a thousand kilometers per second, driven from galax- 
ies experiencing intense recent or ongoing star-formation, 



i.e. starburst galaxies. These superwinds are common in 
the local Universe, occurring in nearl y all galaxies classi- 
fied as undergoing s tarburst activity IlLehnert fc Heckmanl 
ll995llHeckmanll99gT) . and appear ubiquitous among the Ly- 
man break galaxies at redshift ~ 3 where they are blowing 
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~ 100 kpc sized holes in the inter-galac tic medium (e.g. 
IPettini et af]l200ll : lAdelberger et al.ll2003ft . 

Superwinds are believed to be driven by the thermal 
and ram pressure of an initially very hot (T ~ 10 s K), high 
pressure (P/k ~ 10 7 Kcm -3 ) and low density wind, it- 
self created from the merged remnants of very large num- 
bers of core-collapse supernovae (SNe), and to a lesser ex- 
tent the stellar winds from the massive stars, that occur 
over the ~ 100 Myr dur ation of a typical starburst event 
JChevalier fc Clegsl fl985h . The thermalized SN and stellar 
wind ejecta predicted by this model (which we shall call 
the wind fluid for convenience) is too hot and too tenu- 
ous to be easily obs erved (unless heavily mass-loaded, e.g. 
ISuchkov et al] fl996l. but hydrodynamical models of super- 
winds show that the wind fluid sweeps up and incorpo- 
rates larger masses of ambient galactic disk and halo in- 
terstellar medium (ISM) into the superw ind, material which 
is more easily detected observationally ijChevalierfc C leed 
ll985l:ISuchkov et alJll99llStrickland fc Stevensll2000l) . 

Indeed, the majority of nearby superwinds have been 
discovered using optical imaging and spectroscopy to iden- 
tify outflow in the warm ionized gas (gas at T ~ 10 4 K), in 
many cases directly imaging bipolar structures aligned with 
the host galaxy minor axis with kinematics indicative of out- 



Axcji^^T£rvlor]^)78MMcCfu^rv 


1987 


Heckman, Armus, & Milevll987U1990: Bland & Tullv 


198S 


Hcckman. Armus. & Milevll990l:lLehnert & Heckmad 


1995 


19961). Superwinds from galaxies at high redshift are 



recognized by virtue of blue-shifted interstellar absorption 
lines from warm-neutral and warm-ionized gas species (e.; 
IPettini et aTfeood: iFrve et all2002l lAdelberger et alJl2003 
absorption features very similar to those seen in local star- 
burst gala xies with superwind s (Phillips 1993; iKunth et alJ 
119981: iGonzalez Delgado et aDll998t IHeckman et al]l2000D . 
Nearby superwinds have also been extensively studied in 
soft X-ray emission (e.g. jRead et al.l fl997t iDahlem et alJ 



ll998l:IStrickland et all2004ft which probes emission from hot 
gas with temperatures in the range T ~ 10 6 to 10 7 K, fol- 
lowing their initial detection by Ein s tein X-ray observatory 
JWatson. Staneer. fc Griffiths! Il984 iFabbiano fc Trinchieril 
ll^8^ Indeed, observations have demonstrated that all 
phases of the ISM found in normal late type galaxies are 
also incorpora ted into starburst-driven superwinds ( e.g., see 
the reviews of Bland-Hawthorn 1995; Dahlcm 1997). 

Over the last 5 years new observations of superwinds, 
in particular spaced-based observations of Ovi absorp- 
tion and emission from T ~ 10 5,5 K gas in the FUV, 
and thermal emission from T ~ 10 6_7 K gas in X-ray 
regime, have substantially added to our und erstanding of the 
multi- wavelength prop erties o f supe r winds ([ Strkkland e^lJ 
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These observations have shown that the majority of the 
soft X-ray emission in superwind is not due a volume-filling 
wind fluid, but that it and the Ovi-emitting and absorb- 
ing gas arise in some form of interaction between the wind 
fluid and denser ambient disk or halo ISM. This is in qual- 
itative agreement with some of the analytical and theoreti- 
cal models mentioned above. The observations, in particular 
those of Ovi, are not easily explained by the simple ana- 
lytical models of wind-blown bubbles iWeaver et aljlT977t 



iMac Low fc McCrav| ll988l) often applied to superwinds (see 
e.g. IHeckman et aUhoOlbl). 

IChevalier fc Cleeel (I1985T) first posited that both the op- 
tical nebular emission and the soft X-ray emission in su- 
perwinds (as had b een observed by I Axon fc Tavlorlll978l : 
IWatson et alJ 0^984, for example) is from shocks driven 
into clouds embedded in a high velocity wind of merged 
SN ejecta, primarily as the observed soft X-ray luminosities 
exceeded the predicted emission for the wind fluid. Later op- 
tical and low-spatial resolution X-ray observations were in- 
terpreted in the light of interactions between the wind fluid 
and ambient ISM, either with entrained cloud or with am- 
bient disk or halo ISM at the walls of the outflow cavity 
jMcCarthv et al.lll987l : IHeckman et alJll99Ct ISuchkov et"aJJ 
Il994h . In this picture clouds are either dense pre-existing 
structure over-run and incorporated into the wind, or gen- 
erated at the cavity walls through the action of instabil- 
ities in the shell of swept-up ambient ISM (primarily ei- 
ther shell fragmentation through Rayleigh Taylor instabili- 
ties when the superbubble blows out of the disk, or dense 
gas ripped off the walls of the outflow cavity by Kelvin- 
Hclmholtz instabilities). High resolution optical observa- 
tions of superwinds show small scale clumps in the Ha 
filaments of NGC 3079 and M82 with sizes of ~ 20 - 
30 pc, often associated with fainter elongated structures 
JCecil et all200ltlOhvama et alJl2002T) . Thes e are similar to 
the tadpole-like clouds found in s imulations (ISuchkov et alJ 
11994 IStrickland fc Stevensll2000l) . The basic wind/ISM in- 
teraction model for soft X-ray emission in superwinds is 
validated by mo dern high r esolution X -ray observations 
IStrickland et alJbOOCt 120021: IStrickland! I200A ICecil et alJ 
12001 ISchurch et alJ |2002t IStrickland et alJ |2004 . which 
show the soft X-ray emission is structured very similarly 
to optical Ha emission on scales from kpc down to as small 
as 10 - 20 pc, with filaments, clumps and apparently limb- 
brightened walls. Much of this Ha and soft X-ray structure 
is associated with large-scale filaments or the edges of the 
superwind, but examples of wind interactions with clouds 
or larger-scale obstacles exist, e.g. most-notably the M8 2 
northern cloud foevine fc Ballvlll999l: iLehnert et al.lll999l) . 

Returning to the issue of Ovi absorption and emission 
in superwinds, a wind/cloud interaction model is considered 
most plausible for explaining the absorption properties of 
neutr al, warm photoionized and coronal phase gas in NGC 
1705 iHeckman et alJl2001bf) . The kinematics and column 
density of the Ovi-absorbing gas are inconsistent with the 
standard superbubble model, where the coronal phase gas 
arises in a conduction front at the walls of t he superbubble 
shell dCastor et alJll975HWeaver et al.lll977l) . In NGC 1705, 
the superbubble shell (traced by warm photoionized gas) 
has almost certainly ruptured, and is expanding more slowly 
than the Ovi-absorbing gas. The most likely explanation of 
the data is that the Ovi-absorbing gas is being generated as 
the hot gas from the interior of the bubble flows out past 
fragments of the ruptured superbubble shell. 

In this paper we present high-resolution hydrodynami- 
cal models of warm ionized clouds embedded in a superwind, 
and compare the simulated Ovi and soft X-ray properties to 
the existing observational data. The interaction of a super- 
wind with an embedded cool cloud is not the only form of 
wind/ISM interaction considered important in superwinds, 
but wind/cloud interactions undoubtedly do occur within 
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superwinds and we shall focus on them in this paper. These 
models include the effect of thermal conduction, which we 
show plays an important role in shaping both the dynamics 
and radiative properties of the resulting wind/cloud interac- 
tion. Thermal conduction has long been suspected as being 
an important process in superwinds, but has not been in- 
cluded in hydrodynamica l models of superwinds (wi th the 
exception of one model bv lD'Ercole fc Brighentill 999). Sim- 
ulations of cool clouds embedded in a hot fluid, including the 
effects of thermal conduction, have been performed before 
jFerrara fe Shchekinovlll993 IVieser fc Henslerll2000t) . How- 
ever, the models we present here are the first high-resolution 
simulations of an ISM cloud embedded in a hot, supersonic, 
high ram-pressure wind. 

The observational diagnostics available from any sin- 
gle waveband study of superwinds 1 are too crude to 
strongly constrain among the many possible hypotheses 
available. However, multi-wavelength comparisons, e.g. the 
spatial location of Ha emission compared to the soft 
X-ray emiss ion, have proved mor e useful and more ro- 
bust (e.g. iMcCarthv et alj Il987l : iHeckman et alJ Il990t 
IStrickland et, alJl200ft 12002^ . The coronal phase gas respon- 
sible for Ovi emission and absorption in superwinds must ei- 
ther originate from the cooling of originally hotter gas, or is 
collisionally or conductively heated. In any case, Ovi is also 
a tracer of energetic collisional processes in superwi nds in 
the same way X-ray emission is (IHeckman et al]l20021b This 
is unlike the optical nebula emission in superwinds (e.g. Ha 
and [Nil] emission), which can be photo-excited by massive 
stars or from X-ray irradiation, in addition to collisionally 
excited in shocks. The effect of extinction by interstellar ma- 
terial at the wavelengths of the Ovi doublet is equivalent to 
the extinction experienced by soft X-rays of energy ~ 0.3 
keV, somewhat mitigating the observational uncertainties 
in the absolute level of extinction. Oxygen is also the main 
coolant for plasma with temperature between 10 s and 10 6 K, 
and second only in importance to Iron as a coolant in hotter 
gasses. The Ovi-absorbing or emitting gas is also the highest 
temperature phase in superwinds for which the velocity can 
currently be measured, yielding important clues to the rel- 
ative kinematics of di fferent phases in superwinds (see e.g. 
Hec kman et al J2001rj| . Thus a comparison between the Ovi 
and soft X-ray properties of superwinds should hopefully 
prove as, if not more, instructive than the Ha/X-ray com- 
parisons that have already been made. We have thus chosen 
to consider simultaneously Ovi and the soft X-ray obser- 



1 For example, the diagnostic parameters that can be most ro- 
bustly extracted from observation data are: Soft X-ray emission 
— X-ray flux and surface brightness, emission-weighted mean 
temperature, and the apparent gas-phase abundance of O (and 
some of the other alpha-elements) with respect to Fe. Kinemat- 
ics and absolute element abundances, and associated parameters 
such as the emission integral can not currently be obtained or are 
too uncertain to make use of; Ovi emission — fluxes, line-widths 
and simple kinematics; Optical emission — fluxes and surface- 
brightnesses, spatially mapped kinematics and standard nebular 
diagnostics within the inner brighter regions of winds. The nebu- 
lar emission appears dominated by stellar photo-ionization in the 
inner region of winds, but in many cases line ratios (in particular 
Hcr/[N II]) become more "shock-like" at larger distances from the 
central starburst. 



vational properties of our modeled wind/cloud interactions, 
in order to better constrain which models, if any, provide a 
good match to the observed properties of superwinds. 

In Section [5] we describe the numerical code employed 
to perform the models, the parameters of our chosen super- 
wind and cloud model and how we calculate the Ovi and 
X-ray properties thereof. The results of the simulations per- 
formed are described in Section^] The implications of these 
results for understanding superwinds are discussed in Sec- 
tion 0] followed by a summary of the results in Section [£] 



2 THE NUMERICAL METHOD 

As a first step in investigating wind/ISM interactions in su- 
perwinds, we explore the properties of simple model in which 
a single dense cloud (initially at rest) is embedded into a su- 
personic wind of hot gas. 

2.1 Simulations 

To run the simulations we use the 2D BOH (Bologna Hydro- 
dynamics) hydro-code implemented with the thermal con- 
duction. The code is based o n a second-order upwind scheme 
(Bcdosni & D'Ercol elll986l). in w hich consistent advection 
(N orman. Wilson, fc Bartonl l980) is implemented to reduce 
numerical diffusion. We solve the usual hydrodynamic con- 
tinuity equations. We also include a further tracer variable 
representative of the cloud material which is passively ad- 
vected; such a tracer allows us to compute the degree of 
mixing between the cloud and superwind material. 

To take into account the thermal conduction we adopt 
the operator splitting method. We isolate the heat dif- 
fusion term in the energy equation and solve the heat 
transport equation, alternatively along the z and R di- 
rection separately, through the Crank-Nicholson method 
which is unconditionally stable and second order accurate. 
The resulting system of implicit finite difference equations 
is s olved according to the two -stage rec ursion procedure 
(e.g.lRitchmver fc Mortonll967l) . Following rCowie fc McKed 
(1977), we adopt saturated fluxes to avoid unphysical heat 
transport in presence of steep temperature gradients (see 
below). 

In all simulations we adopt 2D cylindrical coordinates. 
The central area of the computational domain is covered by 
an uniform grid with mesh size AT? = Az — 0.1 pc. Beyond 
some distance to the center the linear mesh size increases 
geometrically in both directions with a size ratio of 1.07 
between adjacent zones. The grid edges are thus at large 
distances in order to avoid that possible spurious perturba- 
tions originating at the boundaries may affect the solution 
in the central region of interest. 

For the RE models (see Section [2.21 the total number 
of mesh points is 800 x 400 (z x R), while the uniform grid 
region is covered by of 600 x 300 points. For the models with 
ram pressure the whole grid has 1200 x 400 (z x R) mesh 
points i.e. an extension of ~ 3 kpc in length by ~ 1.3 kpc in 
radius, while the uniform region has 1000 x 300 points. 

In all the models reflecting boundary conditions are en- 
forced on the z-axis. In the RE models, outflow conditions 
are applied at all the remaining boundaries. For the models 
with wind velocity w w 7^ 0, the superwind flows parallel to 
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Table 1. Model parameters for the primary set of superwind/cloud interaction simulations. 



Model 


Pw 




1 w 




X 


M 


Tbs 


CO 


Tec 




(g cm" 


■ 3 ) 


(K) 


(km s _1 ) 






(K) 




(10 5 yr) 


TILP(NC) 


1 X 10- 


-26 


1 x 10 6 


447 


100 


2.98 


2.8 x 10 6 


0.13 


3.3 


TIHP(NC) 


1 x 10" 


-26 


1 x 10 6 


1000 


100 


6.66 


1.4 x 10 7 


3.27 


1.5 


T5LP(NC) 


2 x 10" 


-27 


5 x 10 6 


1000 


500 


2.98 


1.4 x 10 7 


16.34 


3.3 


T5HP(NC) 


2 x 10" 


-27 


5 x 10 6 


2236 


500 


6.66 


6.9 x 10 7 


396.98 


1.5 



Model parameters are the same in models with and without (NC) heat conduction. See Section 12.21 
for details. Parameters with subscript w refer to the undisturbed wind fluid. The density contrast 
between the cloud and the undisturbed wind is x = Pc/fvr- The wind is super-or-hypersonic with 
respect to the cloud, with Mach number M. = i> w /c s ,w The ratio of wind ram pressure to initial cloud 
or wind thermal pressure is ~ 15 in the low pressure (LP) models, and ~ 75 in the high pressure (HP) 
models. The highest temperature in the wind bow shock around each cloud will be approximately 
the post-shock temperature for a strong shock, i.e. Tbs ~ S/lGp/mjiV^/k. The dimensionless heat 
flux saturation parameter <7rj (equation [2J corresponding to this bow shock temperature (assuming 
the bow shock density pes = 4p w ) is a relative measure of the initial strength of the dynamical 
effects seen in the conductive models. The time scale for the cloud to be crushed by the shock driven 
into it is t cc (equation 11*1, 



the z— axis and enters the grid from right-side of the grid 
where inflow boundary conditions are applied. 

The simulations are run until the cloud, which is accel- 
erated by its interaction with the wind, leaves the computa- 
tional grid after ~ 1 to 1.5 Myr. 

Radiative energy losses within the simulation are taken 
into account considering the c ooling curve A for Sola r metal- 
licities as parameterized by iMathews fc Bregmanl lll97Sl) . 
Models of stationary conductive clouds run assuming pri- 
mordial abundances give essentially the same dynamical re- 
sults. With less cooling, there is slightly more hot gas present 
(~ 10% more by mass), but essentially the soft X-ray and 
Ovi luminosity, and Ovi column densities are lower in pro- 
portion to the reduction in metallicity. Clearly in these sim- 
ulations with strong conduction (i.e. at Spitzer levels) the 
structure of conductive front is not controlled or significantly 
influenced by radiative cooling, as in that case we would have 
expected a very significant increase in the total mass of hot 
gas in the low abundance simulations. 

Ovi emission originates in a thin layer of the conduc- 
tive front where the gas is at temperature T ~ 3 x 10 5 K. 
It is crucial to have a numerical resolution good enough to 
capture this thin layer in the conductive front. With a cell 
size of 0.1 pc the Ovi emission from our simulations origi- 
nates from a layer 6-7 grid cells wide. If we reduce the grid 
resolution by a factor of two (to 0.2 pc) cooling rates are 
higher, although the increase in Lovi and Lx is a factor 
of ^ 20% compared to 0.1 pc resolution simulations. This 
suggests that even at 0.2 pc resolution the conductive front 
is still marginally resolved at the Ovi layer, so that we are 
confident that the default 0.1 pc resolution simulations are 
numerical resolving the conductive interface. 

2.2 The wind/cloud model 

We make the following simplifying assumptions: 

(i) The clouds are initially spherical. 

(ii) The clouds are initially in pressure equilibrium with 
the ambient gas. 

(iii) The cloud self-gravity is neglected. 



(iv) Magnetic fields are neglected. 

(v) The gas is not allowed to cool below T = 10 4 K. 

(vi) Non-equilibrium ionization effects on the cooling 
curve are neglected. 

(vii) The elemental abundances of both the wind and 
cloud material are Solar. 

Assumptions (i) and (ii) are made to make simpler ini- 
tial conditions. However, they do not influence very much 
the results because, as it will be described in Section 3, dy- 
namical effects will quickly destroy both the spherical shape 
and the pressure equilibrium of the cloud. About assumption 
(iii), it is easy matter to show that, for the cloud parameters 
adopted in this paper (see below) the cloud mass results al- 
ways lower than the cloud Jeans mass, e ven comput ed taking 
into account external pressure (see e.g.|Shu||l992). In spite 
assumption (iv), we actually consider in Section 3.7 the pos- 
sibility of a reduced value of the thermal conductivity due 
to the presence of magnetic fields. Point (v) assumes that 
the ionizing radiation field due to the massive stars in the 
stellar burst is able to hold the cloud ionized. Assumptions 
(vi) and (vii) are discussed in Section 12.31 and Section 14.31 
respectively. 

The cloud properties, and the range in wind prop- 
erties, we consider are consistent with observational and 
numerical studies of superwinds. At the location of the 
FUSE apert ure used to survey Ovi emission in M82 and 
NGC 3079 jHoopes et all 12003! Hoopes et a l in p repa- 
ration), the Hq-o bservations of lOhvama et all i2002l) and 
ICecil et alJ (l200lT) show ~ 20 - 30 pc-scale clumps within 
the kpc-long filaments. The cloud-like structures visible 
in the highest resolution h ydrodynamical simulations of 
IStrickland fc Stevenj i2000tl have diameters of ~ 30 pc, 
mass density ~ 10 _24 gcm -3 , and are embedded in wind 
fluid with a range of wind velocity, density and ram pres- 
sure similar to those we assume. 

We run different models of an evaporating cloud embed- 
ded into a galactic wind. In most of the models the cloud 
properties are the same: radius R c — 15 pc, temperature 
T c = 10 4 K and cloud mass density p c = 10 _24 gcm -3 
(proton number density n c = 0.42 cm -3 ). The cloud mass 
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Figure 1. (a) Superwind terminal velocity (in km/s) as a function of the efficiency of supernova energy thermalization e and the factor 
by which the hot wind fluid is mass-loaded f3. (b) Oxygen and Iron abundance in units of standard Solar abundance as a function of f), 
assuming the ambient medium has Solar abundances. See Section 12.21 for details. 



is ~ 210Mq. We consider two sets of values for the den- 
sity n w and the temperature T w of the superwind, both 
having the same thermal pressure P — 1.3 x 10 



cm 
10" 



'w, x WJ — v ±u 26 gem 3 
27 g cm" 3 , 5 x 10 6 K). The cloud is thus initially in ther- 



(pw,T w ) = (10- 



10 b K), (p„,T w ) 



12 dyne 
= (2 x 



mal pressure equilibrium with the wind in all the mod- 
els. We further considered two possible values of the ram 
pressure exerted by the superwind on the cloud, namely 
Pram = 2 x 10" 11 dyne cm" 2 (low ram pressure models, 
model name suffix LP), and P ram = 10" 10 dyne cm" 2 (high 
ram pressure models, model name suffix HP). 

Each model is specified by the values of superwind tem- 
perature and ram pressure: for instance, the model T5LP is 
characterized by (n w ,T„) = (8.3 x 10" 4 cm" 3 , 5 x 10 6 K) 
and P r am = 2 x 10" 11 dyne cm" 2 or, equivalently, a super- 
wind velocity with respect to the cloud of « w = 1000 km s" 1 
(see Table 0. To elucidate the role of the cloud size we also 
ran a model equivalent to model T1LP, but with a cloud 
radius R c = 45 pc and with a factor 3 lower resolution. 

For convenience we have chosen to treat the density, 
velocity and metal abundance of the wind fluid as sepa- 
rate variables whose values are not correlated. However, 
theoretically these variables are not independent of each 
other, as mass loaded win ds will have lower v elocity and 
lower metallicity (see e.g. ISuchkov et alJ 1X9961 . To place 
the wind model parameters shown in Table in context 
we show the wind terminal velocity and Oxygen and Iron 
abundances as a function of SN energy thermalization ef- 
ficiency e and mass-loading (3 in Fig. The values shown 
were calculated using the energy and mass return rates and 
elemental yield s from SNe and stellar winds from version 4 
of Starburst99 llLeitherer et al]ll999l) . We assumed contin- 
uous star- formation of IMq / yr at a time 30 Myr after the 
first star formation from gas of Solar abundances, but the 
variation in w w or Z w from assuming a different star forma- 
tion history, time or initial metallicity is typically < 30%. 
The mass injection rate in the starburst M — ft x Msn, 
where Ms at is the mass injection rate due to SNe and stel- 



lar winds alone and /3 is the degree of mass-loading. The 
termin al velocity in the hot wind fluid in lChevalier fc Cleed 
Il985ft superwind model is v x = (2P/M) 5 « 2800(e//3) - 5 . 
To first order the density of the wind fluid scales as /3 1 ' 5 . 
If the elemental abundance of any particular element i is 
Zi.sN in the SN and stellar wind ejecta, and Zijsm in 
the ambient ISM, then metal abundance in the wind fluid 
will be Zj.w = 0' 1 x (Zi.s N + |j9 — 1] X Zijsu)- On 
the lAnders fc Grevessd il989T) abundance scale Zo,sn ~ 
11.5Zo,0 and Z F e,SN ~ 5.2Z Fe:Q . 

In order to disentangle the effect of the heat conduction 
from that due to the ram pressure, we also ran models with- 
out heat conduction (model suffix NC) as well as models 
of evaporating clouds in a static hot medium (model pre- 
fix RE). The model analogous to T5LP, but without heat 
conduction, is called T5LPNC. This primary set of eight 
models are summarized in Table Q We also ran two mod- 
els equivalent to the model with conduction T5LP, but with 
values of the coefficient of thermal conductivity that have 
been reduced by factors 5 and 25 below the Spitzer value. 
A further model equivalent to model T1LP except with a 
cloud radius of R c — 45 pc was performed to investigate 
how the wind/cloud interaction depends on cloud size. The 
static medium models are characterized only by the super- 
wind temperature, thus, for instance, the model analogous 
to T5LP is called RE05 (see Table 0. 



2.3 Analysis 

Calculations of the luminosity, 2-dimensional volume emis- 
sivities, and Ovi column densities from these models we 
performed separately from the simulations themselves. As 
stated above, we assume that the plasma is in collisional 
ionization equilibrium. The electron and ion temperatures in 
each computational cell are thus T = (/x mu P) /(k p), where 
k is the Boltzmann constant, P the thermal pressure, p the 
total mass density, and prriH is the mean mass per particle. 
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Table 2. Model parameters for models of a cold cloud in a sta- 
tionary hot medium. 



Model 


Pw 




1 w 


X 


CO 




(g cm" 


3 ) 


(K) 






RE01 


1 X 10" 


-26 


1 x 10 6 


100 


0.07 


RE05 


2 x 10" 


-27 


5 x 10 6 


500 


8.59 


RE10 


1 x 10" 


-27 


1 x 10 7 


1000 


68.70 



Parameters with subscript w refer to the static hot medium sur- 
rounding the cloud. The density contrast between the cloud and 
the hot medium is x = Pc/pw The dimensionless heat flux satu- 
ration parameter ao is described in Section lij.il equation |21 

For the highly-ionized coronal and hot gas phases we are 
interested in we use a value of fimn = 1.02 x 10~ 24 gem" 3 . 

In the low density limit the power emitted in some en- 
ergy band E, in units of ergs" 1 , from a elemental vol- 
ume SV is 5L E = n e n K A E (T,Z)6V, where A E (T,Z) is 
the emissivity at temperature T, for metal-abundance Z 
in energy band E. For a given mass density p and Solar 
abundances, then to good accuracy % = Xup/mu and 
n e — 0.5 (1 + Xn)p/rnH, where Xu is the Hydrogen mass 
fraction. We use a value of Ah = 0.7057. Making use of 
the assumed symmetry around the z-axis, we sum the lu- 
minosity over all cells within a cylindrical volume chosen 
to encompass the cloud and the majority of the wind bow- 
shock and any cloud fragments. We ignore any absorption 
of this emission from material on the computational grid. 

For the purposes of displaying which regions are the 
strongest Ovi and soft X-ray emitters we plot 2-dimensional 
maps of volume emissivity, i.e. emission per unit volume 
5Le/SV = n e Tin Ae(T, Z), in units of ergs" 1 cm -3 . 

The Ovi emi ssivities we use are based on the MEKAL 
hot plasma code jMewe^^lJh 985f [Ka astra fc Mewelll993r 
iMewe et aflll995l: iLiedahl et alJll99fJl . Note that the lumi- 
nosity quoted is the sum of the two lines in the A = 1032 
A and 1038 A doublet. 

The soft X-ray luminosities we quote are in the E — 0.3 
- 2.0 keV energy band, chosen to correspond to the energy 
band used in the Chandra ACIS observations we compare to. 
The X-ray emissivities used are based on the 1993 update 
to the lRavmond fc Smith! Jl977f) hot plasma code. Although 
the Raymond & Smith code is no longer suitable for detailed 
X-ray spectroscopy, the broad-band emissivities it predicts 
are reasonably accurate. In the chosen energy band the dif- 
ference between the broad-band soft X-ray emissivities given 
by the Raymond & Smith code and the MEKAL code are 
~ 20% for gas with 10 6 < T < 10 7 K, and within a few 
percent for hotter gas. 

Given that the physical size of the observational FUSE 
30" aperture depends on the distance to the observed galax- 
ies (equivalent to ~ 530 pc for M82, and for ~ 1.4 kpc for 
NGC 891, the closest and more distant galaxies we compare 
to), and that we do not know how many clouds actually 
lie within these regions, our quantitative comparison to the 
observational data will concentrate on the ratio of Ovi to 
soft X-ray emission, rather than the absolute emitted lumi- 
nosities. The physical volume represented by the full com- 
putational grid is much larger than the volume we would 
realistically expected to be occupied by only a single cloud 
in a superwind, and the X-ray emission from this volume is 



dominated by the wind, rather than wind/cloud interaction. 
We therefore calculate the emission within a smaller cylin- 
drical region judged large enough to encompass the cloud 
and the most luminous part of the bow shock. This region 
extends 20 pc upstream and 80 pc downstream of the initial 
center of the cloud, and out to a radius of 50 pc. 

We also calculate the predicted Ovi column densities 
through the models, using an upda ted version of the code 
first used in Heckman et al ] feOOlbl) . Along any incremen- 
tal path length Al, the incremental Ovi column density 
AAovi = A/ riH /ovi(T) Ao, where nn is the local hydro- 
gen number density and Ao is the fractional abundance of 
oxygen atoms with respect to hydrogen (Ap = 8.51 x 10~ 4 
for Solar abundance in the lAnders fc Grevesseiri989l scale). 
In collisional ionization equilibrium, and at the low densities 
found in our simulations, the fraction of oxygen atoms in the 
Ovi state is purely a function of temperature .fovi(r), th e 
values for which we take from ISutherland fc Dopital lll993j) . 

The column density is evaluated along between 10 and 
120 lines of sight through each simulation volume. For con- 
venience we calculate lines of sight parallel to the z-axis 
(i.e. parallel to the wind) and within a radius R, which may 
be larger or smaller than the actual radius of the cloud. 
The net Ovi column density quoted is the appropriately- 
weighted average Ovi column density over all lines of sight. 
Column densities calculated using larger radii are typically 
lower than the peak Novi column density that can be found 
using a single line of sight, or even the average over a ra- 
dius R < R c . Obviously the average Ovi column density 
will be reduced if many lines of sight do not intersect any 
gas in the temperature range in which the Ovi ion is abun- 
dant. Within the region of space probed by a typical FUSE 
LWRS aperture there are probably many clouds, but as their 
areal covering factor may in some cases be less than unity 
is worthwhile calculating iVovi over radii R> R c . 

Cooling rates in all the hydrody namical calculations 
assum e Solar abundances (as given in lAnders fc Grevessd 
I1989L which is the commonly used standard in X-ray as- 
tronomy). Although we track the contribution of cloud and 
wind material to the gas density in each computational cell 
and hence calculate the mean metal abundance within each 
cell, for the purpose of calculating the radiative cooling we 
assume Z v = Z c — 1.0 x Zq. In reality cloud and wind 
material are likely to have different metal abundances, with 
the metal abundance of the wind depending on the degree 
of mass-loading. We shall demonstrate later that the rate of 
cooling within the wind is not dynamically significant. This 
allows us to post-facto consider the relative contribution of 
wind and cloud material to the X-ray and Ovi emission in 
cases when Z w 7^ Z c . We shall discuss this issue further in 
Section 0] 

For each computational cell i, j we assess the rel- 
ative contribution by mass from wind material IZij = 



+ Pc 



We assume that the wind and cloud 



material in each cell mix to form a mixture with metallicity 

Zij = (Z w pw,ij 4" Z c p Ci ij^) j \pw,ij 4" Pc.zj) = Z c -\~TLij (Z w Z^). 

Which material dominates the emission from the system is 
then the luminosity- weighted sum of TZij, 



<TL> = 



taking into account the metallicity dependence of the lumi- 
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nosity due to each cell Lij . The mean metal abundance that 
would be observed < Z >, for example in the X-ray band, 
would be the luminosity-weighted sum of Zij. Mathemati- 
cally this is equivalent to < Z >= Z c + < 1Z > (Z„ — Z c ). 

Calculations following the exact ionization balance of 
the wind and cloud material are currently beyond our capa- 
bility. We shall show that radiative cooling is not dynam- 
ically important, but given our interest in the FUV and 
X-ray absorption and emission predicted by these models 
we shall briefly discuss the issue of non-equilibrium effects. 
Non-ionization equilibrium cooling calculations show that 
the emissivity differs most from the equilibrium case for 
gas with a temperature T < 10 6 K. Nevertheless, Ovi col- 
umn densities and Oxygen resonance line emissivities differ 
from the equilibrium values by factors of at most ~ 2-3 
{Edgar fc Chevalieilll986l: ISutherland fc Dopitalll993h . The 
latter paper demonstrates that the Ovi ion is most abundant 
in non-equilibrium cooling conditions at the same kinetic 
temperature as in collision equilibrium. Considering initially 
low-ionization gas being conductively evaporated, the Ovi 
column density m ay be increased by a f actor of ~ 3 from the 
equil ibrium case llWeaver et al.Ml977t llndebetouw fc Shulll 
2004), as the region of Ovi ions extends into lower density 
and kinetically hotter material (the gas is under-ionized, 
in contrast to the over- ionized state in a cooling plasma). 
IWeaver et alJ show that this leads to a minor reduction in 
soft X-ray luminosity, but it is not clear in which direction 
and by how much the Ovi luminosity changes. Thus we be- 
lieve the cooling rates and column densities we calculate 
assuming collisional equilibrium are most likely accurate to 
within a factor ~ 3. This magnitude of effect is not large 
enough to alter our conclusions in a significant manner. 



3 RESULTS 

The main purpose of the present paper is the study of the 
Ovi and X-ray properties of conductively-evaporative cold 
clouds embedded in a hot superwind. However, before to 
present the results obtained for this scenario, it is instructive 
to consider the behaviour of the cloud in two extreme cases: 
i), one in which the evaporating cloud is at rest relative to 
the hot gas, and it), one in which the cloud swept by the 
superwind does not suffer any evaporative mass loss (i.e. 
thermal conduction is completely suppressed). 

3.1 Evaporating clouds at rest in a static hot 
medium (RE models) 

The first approach to the problem of the effects of thermal 
conduction on a coo l cloud embedded in a h ot medium has 
been undertaken by ICowie fc McKed <ll977l) . In the classi- 
cal diffusion approximation of thermal conduction they ob- 
tained an analytical solution for the mass loss rate of an 
evaporating cloud given by 

M cJ = 4.34 x 10- 7 T 6 5/2 i? c , pc M Q yr" 1 , (1) 

where Te, = T/(10 6 K) is the temperature of the ambient gas, 
and i? c ,pc = Rc/(lpc) is the cloud radius. In addition, they 
pointed out that the classical diffusion approximation breaks 
down when the mean free path of electrons becomes larger 
than the temperature scale (saturation). The relevance of 



this effect is quantified by the dimensionless parameter ao 
defined as 

a = 4.22 x 10" 3 — § , (2) 

where nu is the number density of hydrogen atoms in the 
ambient gas (not the number density within the cloud). In 
th e saturated regime (ap > 1) the evaporation rate is given 
bv lCowie et HI ^198 ilk 

M = 1.36A/ C i(j " 5/8 . (3) 

The above results are obtained in the assumption of steady, 
isobaric (and thus low Mach number) flow of the evaporat- 
ing gas, and possible hydrod ynamic effects th a t may arise 
in the system are neglected. IMcKee fc Cowid 1^.975) have 
analyzed the evaporation process in terms of a conduction 
front advancing into the cooler gas. They outlined an in- 
teresting analogy with the ionization fronts, and, just as in 
that case, the velocity of the front v CO nd must be less than 
0.5Cc/c w (D-type front) or greater than 2c w (R-type front), 
where c w is the sound speed of the hot gas, and c c is the 
sound speed of the cloud . In the case of an evaporating cloud, 
ICowie fc McKed (I1977T) have shown that v cond ~ 2aoc^/c w 
for classical evaporation. Thus a critical value <ro,cr = 0.25 
exists beyond which the conduction fronts are D-critical, 
and dynamical effects are important. Saturated evapora- 
tion fronts, instead, have u con d ~ 1.12crJ Cc/c w and are 
always D-critical. In this case the hot pressurized gas behind 
the front expands, driving a shock wave which propagates 
through the cold gas ahead of the conduction front itself. 

The above conclusions apply under the assumption that 
radiative losses are negligible. If such losses exceed the con- 
ductive heat input, material condenses onto the clou d in- 
stead of evaporating. However. IMcKee fc Cowiel (Il977l) have 
shown that this occurs for ao < 0.03, a circumstance never 
met in our models. 

For our models RE01, RE05 and RE10 we have ao = 
0.07, ao = 8.59 and ao = 68.7, respectively. Hydrodynamic 
effects are thus expected to be absent in model RE01, but 
important in models RE05 and RE10. Actually, our simu- 
lations confirm this expectation. In model RE01 the cloud 
evaporates subsonically at the classical rate, without expe- 
riencing any compression. In models RE05 and RE10 shock 
waves develop and compress the evaporating cloud. 

As we show in Section l3.3l all the models of evaporating 
clouds dragged by the superwind develop shock waves driven 
by the conduction front; we thus describe here in detail the 
model RE05, in order to emphasize the cloud dynamics gen- 
erated by the thermal evaporation. As a result of the energy 
transfer by thermal conduction, the pressure at the cloud 
surface rises and drives two opposite flows. The evaporating 
gas moves outward subsonically while the inward flow is su- 
personic (Mach number M > 3) and drives a shock wave. 
This different behavior is tied to the large temperature dif- 
ference between the cloud and the hot gas, and the resulting 
difference in the sound speeds: c w /c c =22.4. 

The spherical shock originating at the surface com- 
presses the cloud as it converges toward the center. The 
evolution of model RE05 is illustrated in Fig. [5] where the 
cloud center is located at r — 0. The first two panels of Fig. [^J 
show the cloud evolution at t = 0.2 Myr and t = 0.4 Myr, 
when the converging shock has not yet reached the center. 
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Figure 2. Density (solid line) and temperature (short dashed 
line) radial profiles as a function of time in the model RE05. 
The cloud center is located at R = 0, while the cloud edge was 
initially at R = 15 pc. The upper two panels show the evolution 
of the spherical shock induced by the thermal conduction while 
it is converging toward the center. The bottom panel displays 
the cloud structure after the shock has bounced back. The long 
dashed line represents the cloud gas density. 



Given the strong radiative losses, the shock is isothermal. 
The cloud temperature remains constant because it is not 
allowed to cool below its initial value of 10 4 K, and thus the 
cloud pressure varies in pace with the cloud density. It is ap- 
parent that the post-shock density (and pressure) increases 
with time, as a consequence of the energy "focusing" as the 
wave approaches the center. This effe ct is expected for a n 
converging adiabatic shock (e.g. lLandau fc Lifshit3ll987f) . 
but remains valid also for a radiative shock, at least in the 
parameters range of our models. 

It is also evident that the density at the outer edge of 
the compressed shell, where the conduction front is located, 
decreases despite the overall compression of the cloud. In 
fact, given the negative value of the radial pressure gradi- 
ent, with time the tail of the density profile lags more and 
more behind; the outer layers of the cloud progressively re- 
tard their collapse, and possibly revert their motion into 
an expansion. The density decrease at the conduction front 
strongly influences the evolution of OVI and X-ray cloud 
luminosities, as discussed in Section f3.1.1l 

At t — 0.43 Myr the shock reaches the center and 
bounces back, re-expanding through the gas still moving to- 




0.5 1 1.5 

t (Myr) 



Figure 3. Time evolution of Lrjvi (solid line) and Lx (dashed 
line) for the RE05 model. The cloud surface area R^ is also shown 
in arbitrary units (long dashed line line). 

ward the center. At this time the central density reaches 
n ~ 2.1 x 10 3 cm" 3 , nearly 5000 times larger than the ini- 
tial cloud density. 

The lowest panel in Fig. shows the density and tem- 
perature profiles at 0.6 Myr. The density jump appearing 
at r ~ 2.5 pc indicates the position of the reflected shock. 
At this stage the central density of the cloud is a factor of 
100 larger than the initial value. The cloud radius attains 
its minimum value R c — 5 pc at t ~ 0.7 Myr, only shortly 
before the reflected shock reaches the edge of the cloud at 
t ~ 0.75 Myr. This latter occurrence influences greatly the 
cloud emission as discussed in the next subsection. 

We point out that, as R c re-expands, the thermal en- 
ergy flux through the cloud surface increases as well. At 
t — 0.95 Myr, when the radius reaches the value of R c = 12 
pc, conditions are restored for a new collapse driven by the 
heat conduction. This oscilla tory behaviour is not found by 
iFerrara fc Shchekinovl l)l993T) who performed similar ID sim- 
ulations; in their models the cloud attains a steady config- 
uration after the initial compression. The difference is due 
to the fact that in their models the cloud is allowed to cool 
to arbitrarily low temperatures, whereas we impose a mini- 
mum temperature to approximate photo-ionization heating 
from massive stars. 

3.1.1 Ovi and X-ray emission in RE models 

As a gas element evaporates from the cloud, it quickly 
goes through a large interval of temperatures, ranging from 
T ~ 3 x 10 5 K where the Ovi emission peaks, up to X-ray 
temperatures (several 10 6 K). Thus significant Ovi and X- 
ray emission originates close to the cloud surface, and both 
are spatially connected. For this reason the time evolution 
of Lovi and Lx are coupled to the cloud dynamics, in par- 
ticular to the variations of the cloud density at R c and to 
the evolution of the cloud size. 

Both Lovi and Lx have similar temporal profiles, both 
coarsely modulated by the behaviour of R c . This is illus- 
trated by Fig. |21 which shows the evolution of Lovi, Lx 
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Figure 4. Time evolution of Lx (upper panel), Lovi (middle 
panel) and log(Lovi/Lx) (lower panel) for all the rest models 
RE. Dashed line: RE01 model; solid line RE05 model; long dashed 
line: RE10 model. 



and L?j? for the model RE05. The sudden initial rise of Lovi 
and Lx is due to the formation of the emitting layer at the 
cloud surface. Successively, the luminosity decreases in pace 
with the cloud radius; we point out, however, that the drop 
of Lovi and Lx can not be explained only on geometrical 
grounds as an effect of the reduction of the emitting surface 
area. Fig. |3] clearly illustrates that the luminosity drop is 
larger than that of R%. Much of the luminosity reduction 
is due to a reduction of the gas density at the conduction 
front, a consequence of the dynamics of the imploding cloud, 
as previously discussed. 

The sudden rise in both X-ray and Ovi luminosity oc- 
curring at t = 0.75 Myr shown in Fig. |3] coincides with the 
arrival of the outwardly- propagating reflected shock at R c . 
This shock compresses the gas of the emitting layer located 
at R c , leading to a rapid increase in emission. 

In general, the evolution of Lovi and Lx is affected 
by a considerable "noise" due to the energy continuously 
injected into the cloud by the conduction front. Even if the 
conditions at the cloud edge are such that a strong shock 
can not be driven by the heat flux, density fluctuations are 
still generated which influence the luminosities. 

The luminosities Lovi and Lx of the three models 
RE01, RE05 and RE10 are compared in Fig.@ In the model 
RE01 the dynamical effects are negligible, and the cloud suf- 
fers a steady evaporation. The rate of such evaporation is 



rather low, M ~ 6.4 x 10 -6 Mq yr _1 , in good agreement to 
the analytical result presented in equation and the cloud 
radius remains constant. The Ovi flux increases slowly with 
time as the conduction front becomes thicker. The X-ray 
flux is not interesting for this model because the tempera- 
ture range through the conduction front is obviously below 
T w = 10 6 , and thus no significant X-ray emission occurs at 
the cloud edge in this case. The constant level of Lx shown 
in Fig. 2] is essentially due to the surrounding hot medium. 

Compared to model RE05, the oscillations of Lovi and 
Lx are accentuated in model RE10, where the heat flux 
is larger because of the larger temperature gradient at the 
cloud edge. 

In the bottom panel of Fig.^Jthe ratio Lovi/Lx is plot- 
ted. We will discuss how this ratio relates to the values ob- 
served in superwinds in Section [I] 

3.2 Non evaporating cloud dragged by the 
superwind (NC models) 

The interaction of a strong shock with a single cloud has 
been the subject of many numerical studies. In the case 
of non-radiative clouds a thorough analysis of this prob- 
lem and a review of relevant literature is provided by 
iKlein et al.l lll994t) . Numerical simula tions in the case of ra- 
diat ive clouds have been presented bv lMellema et ail J2002T) 
and lFragile et all i2004h . In all these papers, after the pas- 
sage of the shock the cloud is embedded in a low-density 
wind, analogous to the situation considered here. However 
the model parameters in the papers just mentioned are not 
appropriate for clouds embedded in superwinds. We thus 
run models of non-evaporating clouds (i.e. without thermal 
conduction) embedded in a superwind. 

The purpose of these models is twofold: on one hand 
they represent an useful check when compared to the results 
by other authors and, on the other hand, they help to un- 
derstand, by comparison with analogous models with heat 
transfer, the role played by the conduction in the general 
dynamics of the cloud. 

In all our models (see Table the superwind is super- 
or-hyper-sonic (Mach number M = 3.2 for models with 
low ram pressure and M = 7.1 for models with high ram 
pressure), despite the high temperature assumed for the 
volume-filling wind material. A bow shock forms around the 
cloud, while a shock is driven into the cloud with velocity 
v s = UwX~ i where \ — Pc/pw The characteristic time r cc 
for the cloud to be crushed by this shock is R c /v s , i.e. 

T CC = X 1/2 i?c/Uw. (4) 

This is the basic timescale governing the evolution of a 
cloud over-run by a blast wave or a wind, wit h cloud de- 
struc tion expected after several crushing times jKlein et alJ 
Il994h . The cloud is fragmented by the action of both Kelvin- 
Helmholtz (K-H) and Ray leigh- Taylor (R-T) instabilities. 
The Richtmyer-Meshkov instability also contributes to the 
fragmentation, but it is less important because it grows lin- 
early rather than exponentially. 

The K-H instabilities are due to the relative motion 
between the cloud and the hot, tenuous background gas. The 
modes destroying the cloud are those whose wavelength is 
comparable to R c , and their growth time IS TKH ~ 'Tec 
Klei n et al J1994I) . This timescale estimate is a lower limit as 
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Figure 5. Density distribution at different times of the cold cloud interacting with the hot tenuous superwind for the reference model 
T5LPNC (no heat conduction). The logarithm of the mass density (in units of gem -3 ) is shown. The distances are given in pc. At the 
beginning of the simulation the cloud center is at z = 20 pc. 



with the formation of the bow shock around the cloud the 
gas actually flowing around the cloud edge has a velocity 
lower than that of the unperturbed superwind. 

The growth time of the R-T instability is trt ~ 
(R c /a) 1/2 , where a ~ irpwR^Vw/Mc is the cloud acceler- 
ation due to the superwind ram pressure. Again, the growth 
time is comparable to the crushing time, trt ~ t cc . Again 
the quoted timescale is a lower limit. In reality the wind 
attempts to form a smooth flow around the cloud, and thus 
the acceleration a is lower than the above estimate. 

The timescales given above hold in the case of an adi- 
abatic evolution of the gas, however they are still ade- 
quate to roughly ch aracterize the flow in the radiative case 
(Frag ile et al.M2004rt . In fact, our radiative simulations sub- 
stantially confirm the picture outlined above. 

In Fig. we show the density evolution of the low ram 
pressure model T5LPNC, for which r cc = 3.3 x 10 5 yr. The 
perturbations grow first close to the axis of symmetry of 
the cloud, where the R-T instabilities are more effective. 
Of course K-H instabilities are also present, and the actual 



growth time is r = (t rt + r KH 



2x1/2 



c/V2 l|Lamblll94it . 



The K-H instabilities are expected to develop faster near 
the tangential point, where the relative velocity is larger. 
However, an inspection to Fig. [£] shows that this region is 
rather smooth; in fact the K-H modes originate at the cloud 
edge but grow downstream, where they are advected by the 
flow. For t < i cc , while the transmitted shock is still crossing 
the cloud, the flow of the superwind results in a low-pressure 
region at the rear of the cloud, into which the cloud material 
moves generating a rarefaction wave which propagates into 
the cloud. This low-pressure region produces also a sort of 



back-flow of the superwind gas, giving rise to a vortex ring 
responsible of the K-H instabilities which severely distort 
the back side of the cloud. 

Although part of the cloud is highly fragmented after 
a few crushing times, its core remains rather compact al- 
though highly distorted. At t = 1 Myr the mass of the core 
is 47 Mq (the original cloud mass was ~ 210A/q), while 
the rest of the cloud, as a consequence of the efficient radia- 
tive cooling, is leaving the grid downstream in form of dense 
filaments. 

Our results are intermediate between those expected 
for a non-radiative cloud and those obtained for a radiative 
cloud. In the former case the cloud would either be com- 
pressed, or be completely destroyed and diffused into the 
ambient medium (Klein et al. 1994). The radiative cloud, in- 
stead, breaks up into numerous, dense, cold fragments which 
survive for many dy namical timescales iMellema et al .120021: 
iFraeile et al.ll2004l) . In our simulations the fragmentation 
process is present but less efficient. This is due to the fact 
that in our models the cloud is not allowed to cool below 
its initial temperature of 10 4 K, while in the papers quoted 
above the minimum allowed temperature is T — 10 K. We 
run a further simulation (not shown here) , analogous to the 
reference model T5LPNC, in which the cloud is now allowed 
to cool down to 10 K. In this case radiative losses are more 
effective and the fragmentation proceeds more efficiently. 
The evoluti on is now very simila r to that found in prev i- 
ous works of lMellema et all l)2002jl and lFragile et all l)2004l) . 
and the cloud is fragmentated in several small clumps after 
few dynamical timescale. 
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Figure 6. Density distribution at different times of the cold cloud interacting with the hot tenuous superwind for the reference model 
T5LP including the effects of heat conduction. The logarithm of the mass density (in units of gem -3 ) is shown. The distances are given 
in pc. At the beginning of the simulation the cloud center is at z = 20 pc. 



3.3 Evaporating clouds dragged by the superwind 

As already discussed, we found that in the RE models 
the dynamical effects linked to the thermal conduction be- 
come important when oo > o- CT . The post-shock gas fac- 
ing the front of the cloud has a temperature Tb s ~ 1-38 x 
10 7 («w/1000 kms -1 ) 2 . As this temperature is always larger 
than the temperature T w of the unperturbed superwind we 
assume, the effective initial value of oo is larger than in the 
static case (see Table0. In particular, <ro is increased above 
(T C i also in the models T1LP and T1HP (which otherwise 
would have uq < <x cr ), and dynamical effects due to the 
thermal conduction are always present in all our models of 
clouds within a wind. 

In Fig. El we show the evolution of the density for the 
reference model T5LP. This model develops a bow shock 
temperature Tb s ~ 1.3 x 10 7 K. As expected, a converging 
shock due to the action of the thermal conduction forms all 
around the cloud (see the shot at t = 0.15 Myr in Fig. 
and the cloud radius has shrunk by a factor of three by 
t — 0.30 Myr. However, the shock formed at the leading 
edge of the cloud is stronger and moves faster through the 
cloud than the shock portion at the rear of the cloud. The 
reason for this is twofold: i) the temperature in front of 
the cloud is larger than at its backside, and ii) the super- 
wind ram pressure contribution is also present at the front 
side. The cloud continues to collapse till ~ 0.4 Myr, when 
the central core reaches a pressure high enough to stop the 
shrinking. At this time the cloud re-expands. The expansion 
occurs preferentially downstream, where it is not contrasted 
by the ram pressure of the incoming superwind. With time 
the cloud looses completely its initial shape and assumes an 



elongated aspect. A successive compression starts at ~ 0.75 
Myr. However, given the complicated cloud morphology and 
the complexity of the flow, the cloud compression is not uni- 
form, but different regions shrink at different rates (see panel 
of Fig. E| at t = 0.8 Myr). After 1 Myr the cloud assumes a 
filamentary shape along the symmetry axis. 

It is worth noting that, contrary to the same model with- 
out conduction, the cloud does not fragment because the R- T 
and K-H instabilities are strongly inhibited by the thermal 
conduction. The K-H instability depends on the the steep- 
ness of the density and velocity gradients across the cloud 
surface. Since heat conduction smoothes o ut these gradients, 
the g rowth of the instabilities is reduced llVieser fc Henslen 
2000). Furthermore the compression of the cloud due to the 
heat conduction reduces the momentum transfer from the 
superwind to the cloud, and thus the cloud acceleration re- 
sponsible of the R-T instabilities. 

In Fig. [7]the evolution of the mass of the cloud and the 
evaporation rate M are plotted for all the models with con- 
duction. This rate is measured by calculating the remaining 
cold gas mass at each simulation time step, specifically the 
mass of gas with T < 1.2 X 10 4 K. In the non-conductive 
models the cloud is always heavily shredded, which makes 
it difficult to measure cloud mass loss, but relatively little 
cold gas gets heated to T > 1.2 x 10 4 K. Focusing on model 
T5LP, at the beginning the value of the mass loss rate is in 
agreement with equation El once the temperature and den- 
sity of the superwind behind the bow-shock is taken into 
account. As expected, the rate drops by a factor of 20 af- 
ter t — 0.35 Myr, consistently with the analogous reduction 
of the evaporating surface area proportional to R^- Then 
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Figure 7. Mass evolution (left panel) and mass loss rate (right panel) for the evaporating cloud models. Black solid lines : T1LP; blue 
short dashed: T1HP; red long dashed: T5LP; green dot dashed: T5HP. 



it increases again following an oscillatory behaviour in pace 
with the evolution of the cloud surface. The reduction of M 
is reflected in the evolution of the cloud mass which shows a 
"flex" in its temporal profile corresponding to the minimum 
of the mass loss rate. The cloud lifetime turns out to be 2- 
3 times longer than that provided by the constant rate M 
given by equation [3] 

The behaviour of all the other models is qualitatively 
similar to that of the reference model. The right panel in 
Fig. |7| reveals that the evaporative mass loss rate is initially 
higher in models with larger values of « w . The same models 
also show a larger frequency in the oscillatory behaviour 
of M. Larger values of t> w lead to larger values of Tbs and 
of the heat flux. The larger the heat flux is, the faster the 
transmitted shock moves into the cloud; as a consequence, 
the re-expansion occurs earlier in models with larger « w . The 
similarity of the mass loss history for models with the same 
w w (at least as long as the clouds maintain a rough spherical 
shape) is rather striking when one compares the temporal 
profiles of M c for models T5LP and T1HP (see the left panel 
of Fig.Q. 



3.4 X-ray and Ovi emission 

3.4-1 Ovi and X-ray emission in NC models 

As pointed out in Section 13.11 when discussing clouds in a 
static hot medium, the interface between the cold dense 
cloud and the hot tenuous wind material is the source of 
both Ovi and soft X-ray emission. In the simulations with- 
out heat conduction these two phases tend to mix at the 
contact surface by numerical diffusion (see Fig. |SJ. While 
such diffusion can mimic qualitatively real phenomena like 
molecular diffusion and turbulent mixing, the quantitative 



effect can be easily overestimated giving rise to excessive 
radiative losses and overly-rapid cloud destruction. 

It has been argued that the numerical spread of the 
contact discontinuities could in principle simulate, at least 
qualitatively, the presence of a heat conduction front; such 
an analogy, however, may be grossly in error in all the cases 
in which the conduction front is expected to give rise to dy- 
namical effects such those described in the previous section, 
and which can not be generated by numerical diffusion. As 
shown in the previous section the inclusion of the thermal 
conduction significantly alters the dynamical behavior of the 
cloud. We also find that the radiative characteristics of the 
cloud change significantly. 

In Fig.|H|we plot the evolution of the Ovi and X-ray vol- 
ume emissivity (i.e. luminosity per unit volume, n e nuA(T)) 
as a function of time for the reference model T5LPNC. In 
the left hand side panels of Fig. [5] the Ovi and X-ray lumi- 
nosities, and X-ray luminosity-weighted wind mass fraction 
< TZ > of all the non-conductive models are compared. 

The Ovi luminosity in model T5LP increases exponen- 
tially during the first ~ 2.4r cc ~ 0.8 Myr (see Fig. [§J, when 
the cloud, although highly distorted, is not yet fragmented. 
The rise in Lovi is due to the increase of the emitting sur- 
face area of the cloud. Subsequently the increase in Lovi is 
even more substantial as cloud fragmentation increases the 
cumulative emitting surface. Finally, after t ~ 1.2 Myr, the 
iovi measured within the simulation volume drops as the 
cloudlets (which are more easily accelerated by the wind) 
start to leave the computational grid. 

Soft X-ray emission originates not only at the cloudlets' 
surfaces but also in the region behind the bow shock, where 
the temperature is T ~ 1.6 x 10 7 K and the density 
~ 2.6 x 10 -3 cm -3 . Wind-based material dominates the X- 
ray emission at times t < 0.8 Myr, based on the X-ray lumi- 
nosity weighted wind mass fraction TZ (bottom left panel of 
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Figure 8. X-ray (left panels) and Ovi (right panels) volume emissivity for the reference model T5LPNC at the same times of Fig. [5] 



The logarithm of the soft X-ray (0.3 - 2.0 keV energy band) and Ovi volume emissivity (in units of ergs 



) is shown. 



Fig. EJ. Once fragmentation becomes significant at t > 0.8 
Myr the fraction of emission from cloud material increases 
to a maximum of ~ 60%. It is worth noting that in Fig. [5] 
the contribution of the X-ray background due to the unper- 
turbed hot superwind is also present. This contribution is of 
the same order of magnitude of that of the cloud, and this 
explains the rather low increase of Lx during the first 0.8 
Myr. Later, analogous to case of Lovi, the X-ray luminosity 
shows a rise and fall due to cloud fragmentation and loss of 
cloudlets from the computational grid. 

Model T1LPNC has the same ram pressure and the 
same analytical crushing timescale r cc of the model T5LPNC 
just discussed, but lower superwind velocity and lower cloud 
to wind density ratio of \ = 100- As pointed out by 
iKlein et alJ lll994) . the time at which the cloud substan- 
tially fragments is also inversely proportional to \- This is 
confirmed by our simulations in which the model T1LPNC 
fragments less dramatically than the model T5LPNC. This 
different dynamical behaviour is responsible for the different 
time evolution of Lovi in model T1LPNC, which increases 
at lower rate and, in particular, does not show the "bump" 
present in the reference model associated to the occurrence 
of the strong fragmentation. Note that the initial value of 
Lovi in model T1LPNC is larger than in model T5LPNC 
because of the different background contribution by the su- 
perwind. The X-ray luminosity in model T1LPNC is always 
dominated by the emission from wind material in the re- 
gion behind the bow shock (1Z is always ~ 1). Here the gas 
has a temperature T ~ 3 x 10 6 K and a number density of 
~ 1.3 x 10~ 2 cm -3 . The density is nearly five times larger 
than in the reference model, and thus Lx is also larger. 

Going back to the high ram pressure models T1HPNC 



and T5HPNC, qualitatively the temporal evolution of Lovi 
and Lx is similar to the low pressure models, although oc- 
curring on shorter timescales. This is as expected as for both 
T1HPNC and T5HPNC r cc = 1.5 x 10 s Myr, shorter than 
in the models with lower ram pressure. The increase in Lovi 
occurs at earlier times because of the faster fragmentation, 
and is higher in magnitude in model T5HPNC because of its 
higher value of x> as discussed above. The strong luminosity 
drop occurring at t — 0.5 Myr is again due to the fact that 
the cloud fragments leave the computational grid. The lumi- 
nosity peak due to the fragmentation is also present in the 
X-ray luminosities, and is associated with a drop in the frac- 
tion of the emission due to wind material. Given the larger 
superwind density in model T1HPNC the bow-shock is of 
relatively greater importance in this model, and the relative 
rise in Lx associated with cloud fragmentation is smaller 
than in T5LPNC or T5HPNC. 



3.4-2 Ovi and X-ray emission of evaporating clouds 
dragged by the superwind 

Fig. 1101 shows the Ovi and X-ray emissivity maps of the 
reference model T5LP at different times. As in the models 
without heat conduction, the X-ray flux originates not only 
close to the cloud surface, but also behind the bow-shock. 
The X-ray and Ovi volume emissivity is largest at the lead- 
ing edge of the cloud, close to the symmetry axis, where the 
density of the evaporating gas has a maximum. This is due 
in part to the larger temperature of the shocked superwind, 
and partially to the compression given by the ram pressure. 

In Fig. |5| (right panels) we plot the time evolution of 
Lx, Lovi, Lqvi/Lx and 1Z for all the models with thermal 
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Figure 9. Time evolution of the logarithm of X-ray (top panel) Ovi (second panel down) luminosities and the logarithm of their ratio 
(third panel down) and X-ray luminosity-weighted wind mass fraction < 1Z > for all models. The right panels refer to models with the 
inclusion of effects due to heat conduction, while left panels refer to non-evaporative models NC. Black solid lines: T1LP; blue short 
dashed: T1HP; red long dashed: T5LP; green dot dashed: T5HP. 
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Figure 10. X-ray (left panels) and Ovi (right panels) volume emissivity for the reference model T5LP at the same times of Fig.[S] The 
logarithm of the soft X-ray (0.3 — 2.0 keV energy band) and Ovi volume emissivity (in units of ergs - 1 cm - ' 3 ) is shown. 



conduction. Concentrating in particular on model T5LP, the 
Ovi and soft X-ray luminosities follow the same oscillating 
behaviour seen in mass loss rates of the cloud (Fig. [TJt , and 
are similar to the analogous model at rest (see e.g. Fig. 
Note that now the luminosities are systematically higher 
than in model RE05, and are more similar to those in model 
RE10. This is due to the larger evaporating flux of emitting 
gas due to the higher bow shock temperatures, as outlined 
above. The low values of 1Z demonstrate that the X-ray emis- 
sion is now dominated by cloud material at essentially all 
times (at least in this case where Z w = Z c ) When compared 
with the T5LPNC model (left hand side panels of Fig. 
both I/x and Lovi are, at least at early times, up to 15- 
25 times higher than the values of the same model without 
conduction. 

The Loyi/Lx ratio in these models is lower than in ei- 
ther the static conductive models, or in the non-conductive 
wind models. This is due to the combined effect of the heat 
conduction and of the bow shock, which tends to increase the 
X-ray flux more than the Ovi flux. It is only in model T1LP 
that conduction does not significantly increase the X-ray 
emission from the wind/cloud interaction. In all other con- 
ductive models X-ray luminosities are higher, and a larger 
and often dominant fraction of the emission comes from ma- 
terial that was originally part of the cloud. 

The time evolution of the luminosities of the other con- 
ductive models is similar to that of the reference model 
T5LP (see Fig. [§J , and the relative differences in the oscil- 
lating period grossly follow the analogous differences found 
in the behaviour of M (see Fig. . In general, the ampli- 
tudes of the luminosity fluctuations are larger for models 
with larger values of Tt, s , and thus larger evaporative M. 



Note, however, that the luminosities in model T1HP reach 
higher values than in model T5LP although this latter de- 
velops essentially the same mass loss in the first 0.5 Myr. 
This is a consequence of the larger value of the ram pres- 
sure which reduces the volume between the cloud and the 
bow shock. This is the region where the volume emissivi- 
ties are largest, and where the evaporating gas in the front 
edge of the cloud moves before to be dragged downstream. 
The increase in the mean density of this region due to the 
reduced volume behind the bow shock thus leads to higher 
luminosities. 



3.5 OVI absorption line column densities and 
kinematics 

In addition to calculating the X-ray and Ovi emission from 
these models we also consider absorption line properties, 
specifically column densities of the Ovi ion probed in FUSE 
obser vations of starbursts (see e.g. iHeckman et alj f2001b. 
I2002D . 

Most of the physical and dynamical effects that alter the 
previously-discussed X-ray and Ovi emission also affect the 
gas that would be seen with absorption line probes. We have 
therefore kept this discussion short, and present only the 
time-averaged Ovi column densities for each of the models 
in Table [I] For each model we report the mean value of 
iVovi through the center of the cloud averaged over two 
different radii: 7? sor = 5 pc and i? S or = 15 pc (this latter 
being equal to the cloud radius). Note that these radii give 
rise to a spatial area lower than that subtended by a FUSE 
LWRS aperture. 

Thus, if the areal covering fraction of clouds within an 
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Table 3. Mean fractions of wind mechanical energy radiated and Lovi/i'X ratios in the wind/cloud interaction, and in the 
undisturbed wind. 



Model 




< Lx/Imec > 


< ioVl/^mcc > 




< Lqvi/Lx > 


^avg 




erg a 


(96) 


(%) 


(%) 


{Z v = Zq) 


(Z w = 10Z Q ) 


Myr 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


TILP(NC) 


1.34 x 10 36 


0.23 (0.122) 


0.20 (0.13) 


0.89 (0.55) 


0.87 (1.07) 


0.1 (0.11) 


1.5 (1.5) 






0.045 


0.028 


0.138 


0.62 


0.62 




TIHP(NC) 


1.50 x 10 37 


0.07 (0.002) 


0.05 (0.11) 


0.25 (0.37) 


0.71 (44.5) 


0.12 (5.23) 


1.0 (0.8) 






0.0040 


0.0025 


0.0123 


0.62 


0.62 




T5LP(NC) 


3.00 x 10 3e 


0.12 (0.013) 


0.13 (0.14) 


0.55 (0.49) 


1.11 (11.2) 


0.34 (1.68) 


1.0 (1.5) 






0.0062 


9.7 x 10~ 7 


0.0062 


1.6 x 10" 4 


1.6 x 10~ 4 




T5HP(NC) 


3.33 x 10 37 


0.03 (0.004) 


0.07 (0.08) 


0.29 (0.26) 


2.68 (18.5) 


1.31 (7.58) 


0.55 (0.8) 






0.00056 


8.8 x 10" 8 


0.00056 


1.6 x 10" 4 


1.6 x 10" 4 





Notes — The time-averaged mean X-ray, Ovi and total luminosities are given as a fraction of the superwind mechanical power 
Lmoc = 0.5 p w v^f 7r -R^m passing through the chosen volume (a cylinder of radius R = 10 pc and extending 20 pc upstream and 80 
pc downstream of the initial cloud center). Note that the luminosities are calculated over a larger volume, with a radius of 50 pc. 
See Section l4.4l for details. Columns are as follows: Column 1: Model name. Values given in parentheses correspond the cloud models 
without thermal conduction. The second line of values associated with each model shows the Lx/L m ec, iovi/^mec, itot/imec 
and -Lovi/^x values for an equal volume of wind material with no clouds within it. Column 2: Wind mechanical luminosity L me c- 
Column 3: Average soft X-ray emission (E=0. 3-2.0 keV) to wind mechanical power ratio (in percent). All values in columns 3 - 
6 were calculated using Z w = Z c = Zq. Column 4: Average Ovi doublet (A = 1032 A and 1038 A) emission to wind mechanical 
power ratio (in percent). Column 5: Average total FUV and soft X-ray em ission to wind mec hanical power ratio (in percent), 
assuming the Ovi emission accounts for 30% off the total FUV cooling rate iHoopes et aljl2003h . Column 6: Average Ovi to soft 
X-ray luminosity ratio. Column 7: As column 6, except calculated with Z c = Zq and Z w = IOZq. Column 8: The time over which 
the luminosities were averaged. 



LWRS is of order unity, then this is the column density that 
would be observed. But, if the cloud covering factor is less, 
then the column density that would be observed with the 
LWRS on FUSE would be proportionally less. 

In the models without conduction the Ovi column den- 
sity A r ovi varies with time, depending on the fragmentation 
of the cloud, and roughly proportional to the behaviour of 
Lovi- From Table 0] it is clear that the time-averaged mean 
column density < iVovi > is higher for models with larger 
wind ram pressure (and consequently shorter cloud crushing 
time Tec). Models with the same ram pressure have roughly 
the same value of iVovi (within 0.1 dex). Analogously to 
our earlier discussion on the Lovi behaviour, the model 
T5LPNC has a slightly higher mean iVovr value compared 
to equivalent ram pressure model T1LPNC, most probably 
because of the higher cloud fragmentation (which depends 
on the cloud density contrast \). An effect of similar mag- 
nitude is seen comparing the non-conductive high pressure 
models. 

These effects are most clearly seen along lines of sight 
weighted toward the center of the original cloud where the 
density of T ~ 10 55 K gas is greatest, i.e. R SOI = 5 pc in 
Table|l] Mean density columns over a larger radius still show 
the same general pattern in these non-conductive models, as 
cool clouds fragment and their Ovi interfaces extend to radii 
greater than that of the original cloud. 

The inclusion of thermal conduction alters the mean 
Ovi columns as follows. Conductive effects crush the cloud, 
and cloud fragmentation is almost completely suppressed. 
As a result there is a large difference between the mean 
iVovi column density averaged over a line of sight of 5 pc 
radius (which always intersects the cloud) to that averaged 
over 15 pc radius (which often includes line of sight with 
negligible Ovi-absorbing material). 

Considering R SOI — 5 pc, the < iVovi > value observed 



is physically driven by the bow shock temperature Tb s (as is 
the case for Lovi), with higher T^ s driving more conductive 
evaporation and hence creating more T ~ 10 5 ' 5 K gas. Note 
how much closer the iVovi values are between models T1HP 
and T5LP (which have very similar post-shock temperatures 
Tbs) than in the case of the equivalent non-conductive mod- 
els. The higher iVovi value for T5LP is due to the higher 
initial temperature of the wind in T5LP. This hotter gas 
also drives stronger evaporation behind the cloud (on the 
reverse side to the wind bow shock) than in model T1HP. 

When averaging the Ovi column over the larger radius 
of -R S or = 15 pc, the time-averaged Novi is relatively uni- 
form, even slightly decreasing as T^ B increases. Here the 
strong increase in Novi with increased 7b s along lines of 
sight passing through the cloud is offset by increased num- 
bers of lines of sight with negligible Nam, as increased Tbs 
leads to stronger cloud compression and reduces its cross- 
section. 

Within the short period of time over which we can fol- 
low the evolution of the cloud in these simulations the cold 
gas is accelerated by the wind to a mass-weighted average 
velocity between 20 and 100 km s _1 (at t = 0.5 Myr). Cold 
gas velocities are systematically larger in the high ram pres- 
sure models compared to low pressure models. Comparing 
equivalent conductive and non-conductive models, the mean 
cold gas velocity is higher in the non-conductive models as 
cloud fragments are easily accelerated to higher velocity. The 
fraction of the wind mechanical energy converted to cloud 
kinetic energy by this time (calculated within a radius of 
10 pc for all models) ranges between 0.6 and 6%, and on 
average is 2% for models with conduction and 5% for non- 
conductive models. Note that more wind mechanical energy 
is converted to cloud kinetic energy than is radiated away 
(see Table EJ. 

We investigated the Ovi kinematics of a randomly- 
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chosen subsample of the simulated Ovi absorption lines. 
Mean velocities were in the range to 200 km s -1 , although 
we did not find any clear correlation between mean cold gas 
velocities and Ovi profile velocities. We found line widths 
equivalent to b values in the range 30 to 200 km s - . 

3.6 The influence of the cloud radius 

In all the previously described simulations, the initial radius 
of the cool cloud was 15 pc, consistent with the smallest 
scale structures observed in the warm ionized gas in super- 
winds and seen in numerical simulations of superwinds (as 
discussed in Section 0. This does not mean that cool clouds 
of only this size exist. Much larger structures of relatively 
cool T < , 10 4 K gas exist (e.g . the northern cloud or cap in 
M82, see lLehnert et al . 1999, or the large scale filaments in 
the halo of NGC 253) although how coherent these struc- 
tures are is hard to assess, and these very large structures 
may actually mark to edges of the flow rather than being 
fully embedded within the wind. 

The distribution of cloud sizes and masses in winds has 
not been determined observationally or theoretically. As the 
average cloud size may differ from the R c = 15 pc used as 
a default in our simulations it is worth considering how the 
emission and absorption properties of a sup erwind /cloud in- 
teraction depend on cloud radius. To investigate this we ran 
one further model, with a cloud of radius _R C = 45 pc (three 
times larger than in all the previously discussed models), on 
a grid with resolution Ar — Az — 0.3 pc. All other model 
parameters were the same as in the low pressure conductive 
model T1LP. 

We compared the soft X-ray luminosity, Ovi luminosity 
and time-averaged Ovi column density in this large radius 
model to the values found in model T1LP, after subtracting 
the background due to the undisturbed superwind within 
the same volume. As described in Section |5] for the model 
with _R C = 15 pc the soft X-ray and Ovi emission was cal- 
culated within a radius of 50 pc off the cloud axis, and 20 
pc upstream and 80 pc downstream of the initial location of 
the cloud. As this is similar to the size of the cloud in the 
larger cloud radius model we also calculated the emission 
from the large cloud radius model within a region 3 times 
larger in each dimension than that used in the default mod- 
els. Column densities were calculated within the same 5 and 
15 pc radii, but along a path length 3 times larger. As the 
size of regions used differ between the default model and the 
large cloud radius model, we subtracted the emission or ab- 
sorption contribution expected from the undisturbed wind 
within the region to allow a fair, if approximate, comparison 
of the relative change in luminosity or column density. 

The soft X-ray emission in the large cloud radius model 
is ~ 40 times greater than in model T1LP, i.e. scaling 
roughly as Lx oc J?|?' 3 . As the soft X-ray emission in this 
supersonic wind/cloud interaction mainly arises in the bow 
shock (in model T1LP conduction plays less of a role in 
enhancing the X-ray emission than in the other conductive 
model), this scaling can be understood as cross section of 
the bow shock scaling as Rj: and its thickness scaling as Rl, 
and thus the luminosity scales as the approximate increase 
in shocked gas volume. 

In contrast to the X-ray emission, which is spread over 
a large volume within the bow shock, the Ovi emission is 
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Figure 11. Time evolution of the logarithm of the X-ray (upper 
panel), Ovi (middle panel) and the logarithm of their ratio (bot- 
tom panel) for model T5LP for different values of the reduction 
factor / of the Spitzer heat conduction coefficient. Balck solid 
line: / = 1; blue short dashed: / = 5; red long dashed: / = 25; 
green dot dashed: / = co. 



concentrated onto a relatively thin interface region at the 
immediate edge of the cloud core. The relative increase in 
Lovi is only a factor ~ 12, which indicates that Lovi oc 
R~ 2 . This suggests that the thickness of Ovi zone of the 
conductive front is relatively independent of the cloud radius 
(to within ~ 30%), and that most of the Ovi luminosity 
increase is due to increased cloud surface area. 

This conclusion is supported by the simulated Ovi col- 
umn densities, which are very similar between the large 
cloud radius model and model T1LP. For the large cloud 
radius model the time averaged mean is log iVovi = 12.89 
for 7? sor = 5 pc (higher by ~ 7%), and logiVovi = 13.11 
for 7? sor = 15 pc (higher by ~ 70%). As in model T1LP the 
cloud is compressed to radii less than 15 pc so that some 
lines of sight see little or no Ovi-absorbing gas, the fairest 
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Table 4. Mean Ovi column densities through the simulated clouds. 



Model 


log 


(JVovi) 


Model 


log 


(AW 




-Rsor = 5 pC 


Rsor = 15 pC 




-Rsor = 5 pC 


-Rsor = 15 pC 


T1LP 


12.86 


12.88 


T1LPNC 


13.08 


12.54 


T1HP 


13.29 


12.83 


T1HPNC 


13.29 


12.98 


T5LP 


13.43 


12.80 


T5LPNC 


13.19 


12.97 


T5HP 


13.52 


12.80 


T5HPNC 


13.38 


13.05 



Notes — Time-averaged mean values of log./Vovi ( m units of cm~ 2 ) due to the 
wind/cloud interaction in each of the models. See Section l3.5l for details. The Ovi col- 
umn density due to the undisturbed wind within the computational volume has been 
evaluated and subtracted from the total derived Ovi column density, i.e. the quoted 
column density is the additional Ovi column generated by the wind/cloud interaction. 
The undisturbed wind has log./Vovi = 12.65 in all Tl models, and log./Vovi = 10.04 
in all T5 models. 



comparison is between the columns averaged over the R = 5 
synthetic beam. The similarity in column density is thus 
consistent with a conductive interface whose density and 
thickness is independent of the cloud radius to first order. 

Thus these results indicate that important emission 
properties, specifically the soft X-ray emission, but also the 
related Lovi/Lx ratio and the fraction of wind energy lost 
to radiation (Lx/L mcc and L to t/L mcc ), do depend on the 
cloud radius (larger clouds will give lower Lovi/Lx ratios 
and higher Lx/L mcc ratios). 



3.7 Lower values of the coefficient of thermal 
conductivity 



Table 5. Mean Ovi to X-ray flux ratios and Ovi column den- 
sities for different values / of the reduction of the coefficient of 
thermal conductivity, in models based on T5LP conditions. All 
calculations assume Z w = Zq. 
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< Lovi/Lyi > 




log (iVovi) 






Rsor — 


5 pC -Rsor = 15 pc 


1 


1.1 


13.43 


12.80 


5 


1.0 


12.39 


12.04 


25 


2.3 


12.11 


12.20 


oo 


11.2 


13.19 


12.97 



Notes — see Section l3.7l for details. 
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2004), suggest that the actual heat conduction coefficient 
is reduced in the presence of tangled magnetic fields by 
a factor of at least 5 relative to the Spitzer value. In 
order to understand how sensitive our results are to the 
absolute value of the conductivity coefficient « we ran two 
more models analogous to model T5LP, but reducing this 
coefficient by factors / = 5 and / = 25 respectively, i.e. 
« = «s P //, where K Sp = 6.1 x 10" 7 T 5/2 erg s K _1 cm -1 
is the "standard" Spitzer value (adopted in Eq. 1 and in 
all the simulations with heat conduction discussed so far. 
Thus / = 1 represents model T5LP, and / = oo is model 
T5LPNC). We also ran a model with / = 125, but the 
results were very similar to those with / = oo, so we do not 
discuss this model further. 

As expected the apparent strength of the R-T and K-H 
instabilities increase as / increases. However, within t = 1 
Myr (the time at which the simulations have been stopped) 
the cloud does not undergo any significant fragmentation 
even for / = 25. This fact influences the temporal behaviour 
of Lx and Lovi- The upper panel of Fig. 1111 illustrates Lx 
temporal profiles of the model T5LP for the four values of 
/, / = (1, 5, 25, oo). The profile for / = 5 is similar in shape 
to that for / = 1 but at lower X-ray luminosity because 
of the weaker conductive mass loading of the bow shock 
region. The oscillations are weaker because the reduced heat 
conduction produces smaller oscillations of the cloud radius. 



The X-ray luminosity for / = 25 is rather constant and only 
slightly larger than that for f — oo. Note however that, for 
t > 0.8 Myr, Lx for the model T5LPNC increases because of 
the fragmentation, while remains nearly constant for / = 25. 
For the same reason, Lovi remains quite low for f — 5 and 
/ = 25 relative to / = oo (second panel of Fig, lilt . 

The cloud lifetime for these models with conduction of 
intermediate strength (5 < / < 25) is increased because 
of the reduced evaporation rate while the fragmentation re- 
mains essentially inhibited. 

In Table |S] we summarize the time- averaged mean ratio 
Lovi/Lx for different values of /. It is interesting to note 
that this value is essentially the same for values of / = 1 
and f — 5, and only increases to 2.3 for / = 25. As shown 
m Fig. [TT] the increase in Lovi/Lx ratios with increasing 
/ is due to the gradual reduction in Lx while Lovi values 
are relatively similar irrespective of the strength of thermal 
conduction. 

It is also interesting to compare the different values 
of the mean OVI column density for different values of / 
(shown in Tabled, where two values of 7? sor are considered. 
For Rs OT — 15 the column density is reduced by a factor of 
5 as / increases from / = 1 to / = 5. This can be easily 
understood because the amount of the OVI emitting gas is 
roughly proportional to the rate of mass loss, and thus to 
/ . We point out that for R SOI = 5 the difference is larger 
because, contrary to the case / = 1, with / = 5 the cloud 
shrinks less and does not fit entirely inside -Rsor- As / in- 
creases further, the column density reverts this trend and 



Conductive clouds in supernova- driven super-winds 19 



increases again. This is due to the effect of the numerical 
diffusion which starts to prevail on the thermal conduction 
in spreading the cloud edge when 5 < / < 25. The column 
density increase for / > 25 is due to the occurrence of the 
fragmentation and the subsequent increase of the emitting 
surface. 

These results suggest that even with thermal conduc- 
tion inhibited by realistic amounts the interesting dynamical 
and radiative effects of conduction on wind/cloud interac- 
tions are retained. X-ray emission is enhanced, the Lovi/^x 
ratio remains low, and hydrodynamical cloud fragmentation 
is inhibited. 



4 DISCUSSION AND IMPLICATIONS FOR 
SUPERWINDS 

4.1 Cloud survival in superwinds 

In order to produce the broad absorption line profiles of 
warm, neutral and coronal gas observed in superwinds, am- 
bient ISM material must survive being entrained within 
the hotter superwind material long enough for it to be ac- 
celerated from rest to veloci t ies of several h undred km/s 
l|Heckman et al J 1200(1 l2001bl iHeckmadEoOsh . Warm ion- 
ized gas is observed to extend out to heights of iifwiM > 
10 kpc from the host galaxy in many superwinds. If this 
material has been transported from an original location 
within the disk this would imply a life time of t ~ 
2 (_HwiM/kpc)/(i>wiM/500 kms -1 ) Myr. Yet clouds over- 
run by a blast wave or enveloped in a supersonic wind are 
expected to be hydrody namically destroy ed within a few 
crushing time scales (see lKlein et al.lll994h . which are typi- 
cally less than 1 Myr for typical estimates of superwind and 
cloud parameters (see Table 0. 

Understanding the conditions most conducive to cloud 
survival is therefore important. Theoretical models in which 
clouds are rapidly destroyed are also less likely to be good 
representations of the physics of superwinds than those in 
which clouds have longer lifetimes. 

In this work we find that cloud models with thermal 
conduction have very different properties in terms of cloud 
survival to those without conduction. Thermal conduction 
suppresses the instabilities that strongly fragment the cloud 
in the non-conductive case, but also leads to significant cloud 
mass loss by evaporation into the hotter phases. 

The mass and mass loss rate in cold gas (log T < 4.08) 
is shown as a function of time in Fig. Q f° r the models 
with thermal conduction. The majority of this mass loss 
is due to conductive-driven evaporation, as cold cloud gas 
is heated to T > 10 5 K. As previously discussed in Section 
13.31 the net evaporation rate is proportional to some power 
of the temperatures of the bow shock around the cloud, so 
the evaporation rate is least in the low wind velocity, low 
wind temperature model T1LP. The saturation of thermal 
conduction reduces the evaporation rates compared to the 
classical case (Equ.|3J, but this analytical theory does not 
include dynamical effects (such as cloud compression) that 
further reduce the evaporation rate in models T1HP, T5LP 
and T5HP. The net effect is for an increase cloud survival 
times in cases where op > 0.25 in com parison to the values 
predicted from lCowie fc McKed lll977ft . 



If we define a cloud lifetime as Tnf e = M c i ou d/M c ioud, 
and using the mean cloud mass loss rate found in the sim- 
ulations, then the cloud life times in the conductive models 
are as follows. Model T1LP: 12.8 Myr; Model T1HP: 3.3 
Myr; Model T5LP: 2.0 Myr; Model T5HP: 0.9 Myr. With 
the exception of model T1LP, these lifetimes are somewhat 
less than the lifetime we expect the average cloud in a su- 
perwind should have. 

The mean mass loss rate in the simulation of the cloud 
with Rc = 45 pc is M c i oud = 6.26 x 10~ 5 Af©yr _1 , three 
times larger than in the analogous model T1LP with the 
standard radius Rc = 15 pc (1.93 x 1O" 5 M yr" 1 ). Con- 
duction is not saturated in this model, as ao = 0.043, so 
that the mass loss rate is classical and scales in propor- 
tion R c . The numerical result is consistent with the the- 
oretical expectation. In fact, the ICowie fc McKeel formula 
(Equ.0 yields M c ioud = 1.84 x l(T 5 Af s yr _1 for R c = 15 
pc, and 5.52 x lO" 5 M0yr _1 for R c = 45 pc (These values 
are obtained considering a temperature T = 2 x 10 6 K mea- 
sured from the simulations between the bow shock and the 
cloud [instead of the theoretical bow-shock temperature], 
and halving the values obtained considering the lower tem- 
perature at the rear of the cloud). The cloud lifetime scales 
as R 2 in the classical regime, as expected. Note however that 
if dynamical effects were important (as in all models but 
T1LP) we expect a lifetime longer than the theoretical one. 
To summarize, cloud lifetime is a strong function of cloud 
size. Conductive clouds larger than the default R c = 15 pc 
we consider may well have survival times equivalent to wind 
flow times, even in the case of conduction at Spitzer levels. 

It is difficult to assess cloud lifetimes in the non- 
conductive models. Cloud fragments are very rapidly accel- 
erated and advected off the fixed numerical grid we have 
used in these simulations. It is hard to assess what fraction 
of the mass loss from the cool phase is due to this process, as 
opposed to transfer of mass into hotter phases, or to algorith- 
mically determine the mass of the remnant of the original 
cloud. Bearing these complications in mind, we have sim- 
ply totaled the mass of all cool gas (again logT < 4.08) on 
the computational grid. For each model, comparison to the 
equivalent conductive model at the same time shows that for 
all cases there is more cool gas in the non-conductive models 
by a factor of typically < 15% in the low wind temperature 
models (Tl models), and ~ 50 to 100% in the high wind 
temperature models (T5 models). Note that accounting for 
cloud fragments lost from the grid would increase these dif- 
ferences. In this crude sense mass loss (to the hot phase) 
is less in the non-conductive models. The higher total cool 
gas mass in the NC models is spread over multiple cloud 
fragments. 

For model T5LPNC we did identify the cloud core by 
hand at t = 1 Myr and obtained a cold gas mass for this 
core of 47M© . At this time the cold gas mass in conductive 
model T5LP is ~ IOOMq. Thus in terms of maximizing the 
lifetime of the original cloud non-conductive models are not 
better than the conductive models. 

The simulations we have presented should be consid- 
ered to as two extremes in which cloud lifetimes are min- 
imized: clouds experiencing conduction-driven evaporation 
due to conduction at effectively saturated levels, verses ef- 
ficient hydrodynamical fragmentation under conditions of 
no conduction at all. Our brief investigation of models with 
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Table 6. Observed FUSE Ovi to soft X-ray flux ratios in superwinds and normal spiral halos. 



Galaxy and region 


Width 


fovi 


fx, 0.3-2.0 


log L x 


/ovi//x 


Ovi ref. 


X-ray ref. 




pc 


ergs -1 cm -2 


ergs -1 cm -2 










(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


M82-A 


525 


< 11.4 x 10 -14 


(94.6 ± 2.2) x 10 -14 


39.2 


< 0.12 


1 


1 


M82-B 


525 


< 11.4 x 10 -14 


(69.5 ± 2.1) x 10 -14 


39.0 


< 0.16 


1 


1 


M82-C 


525 


< 9.3 x 10 -14 


(11.9 ±0.8) x 10 -14 


38.3 


< 0.78 


1 


1 


M82-D 


525 


< 9.5 x 10 -14 


(6.3 ± 0.5) x 10 -14 


38.0 


< 1.51 


1 


1 


NGC 3079 


2490 


< 1.5 x 10 -14 


(30.7 ±2.9) x 10 -14 


40.0 


< 0.05 


2 


2 


NGC 4631-A 


1092 


(3.6 ±0.8) x 10 -15 


(3.6 ± 0.9) x 10 -15 


37.4 


1.00 ±0.33 


3 


4 


NGC 4631-B 


1092 


(6.2 ± 1.0) x 10 -15 


(10.0 ± 1.2) x 10 -15 


37.8 


0.62 ±0.12 


3 


4 


NGC 891-2 


1395 


(2.1 ± 1.1) x 10 -15 


(7.8 ± 3.8) x 10 -15 


37.9 


0.27 ±0.19 


3 


4 



A compilation of observed FUSE Ovi flux measurements in the halos of nearby starburst and normal spiral galaxies, along with 
Chandra ACIS-S measurements of the diffuse soft X-ray luminosity within the same region. The physical size corresponding 
the 30" angular width of the square FUSE LWRS aperture is given in column 2, assuming distances of 3.6, 17.1, 7.5 and 9.6 
Mpc to M82, NGC 3079, NGC 4631 and NGC 891 respectively. The Ovi fluxes quoted in column 3 refer to the sum of the 1032 
and 1038 A lines, assuming the flux in the 1038 A line is exactly half that of the 1032 A fluxes measured (the reference for 
which is given in column 7). The diffuse soft X-ray (E = 0.3 — 2.0 keV) fluxes are based on Chandra ACIS-S measurements, the 
reference for which is given in column 8. The log of the absorption-corrected soft X-ray luminosity (in units of ergs -1 ) within 
the FUSE observation region is given in column 5. Upper limits are 3a. Errors in fluxes, where given, are purely the statistical 
uncertainties in the number of Ovi 1032 A or diffuse E = 0.3 — 2.0 keV X-ray photons within the re gion encompassed by the 
FUSE apert ure. All fluxes are corrected for extinction/absorption, using the Galactic extinction law JCardelli et al.lll98flf) for 
Ovi and the lMorrison fc McCammorJ ll983l) model for X-ray absorption cross sections. Only Galactic foreground extincti on is 
assumed for NGC 4631 and NGC 891, while intrinsic absorption within the wind is required for M82 iHoopes et aljl2003F> and 
NGC 3079 (Hoopes et al, in preparation). 

References: (T) lHoopes et ai1l2003tV The X-ray fluxes are based on spectral fits to the diffuse X-ray emi ssion within each of the 
regions covered by the FUSE apertures. (2) From Hoopes et al (in preparation). (3jJO l tteje^^lJj20Q3j^ 1 _(4) This work. X-ray 
fluxes are based on scaling the measured halo- region diffuse soft X-ray fluxes from IStrick land et al] 120041) by the ratio of the 
E = 0.3 — 2.0 keV diffuse count rate in the region of the FUSE aperture to the total halo diffuse count rate. 



conduction at sub-Spitzer levels shows that in certain cir- 
cumstances conduction is strong enough to suppress the hy- 
drodynamical instabilities that fragment clouds, while also 
reducing cloud evaporation. As previously discussed larger 
clouds also have longer lifetimes. Thus it seems plausible 
that it will be possible to find conditions in which simulated 
cloud lifetimes are the > 10 Myr timescales we expected 
from observations, but further simulations will be required 
to test this hypothesis. 



can wrap around a cloud) than shock crushing, and achieves 
high compression ratios. 

The clouds we have simulated are not natural candi- 
dates for star formation, given their low initial densities and 
masses. Nevertheless, conduction-driven implosion is a very 
interesting mechanism than deserves further investigation, 
and should be considered when thinking of extra-planar star 
formation in actively star-forming galaxies, and in galaxies 
over run by superwinds or jets. 



4.2 Conduction-driven cloud implosion and star 
formation 

The possible role of interstellar or inter-galactic shock waves 
in promoting star formation (or even galaxy formation) 
by crushing am bient gas clouds has been considered many 
times (e.g. see IWoodwardl Il976t lOstriker fc Cowiel 1981; 
Shchekinovl Il989t iBeeelman fc Cioffil Il989l: iBicknell et aD 
20001: iGorti fc Hollenbachl 120021: iFraeile et al.1 l2004 
Sc^am£U3ie^c^^dTl20o3) . Nevertheless, it has long been 
known that shock-crushing does have some undesirable 
effects, as the co mpression is not isotropic but tends to 
pancake the cloud dWoodwardlll976J) . 

The conduction-driven cloud compression seen in 
our si mulations, and first noted by iFerrara fc Shchekinovl 
(1993), should also be active in many of the contexts in 
which shock-driven cloud-crushing is considered, as long as 
hot gas is generated and thermal conduction is not com- 
pletely suppressed. This conduction-driven compression is 
much more isotropic (due to the rapidity with which hot gas 



4.3 The origin and metallicity of the X-ray 
emitting gas 

What implications do these simulations have for understand- 
ing soft X-ray emission in superwinds? As discussed in Sec- 
tion from observations it is clear that the majority of 
the soft X-ray emitting gas comes from a relatively small 
fraction of the total volume occ upied by the wind. Cu rrent 
X-ray observational studies (see IStrickland et al.ll2004T . and 
references therein) are consistent with the following scenar- 
ios: that the hot gas is ambient gas heated by the wind 
to X-ray-emitting temperatures, or alternatively is denser- 
than-average regions of initially hot merged SN and stellar 
wind ejecta, or is some mixture of the two. This is an im- 
portant issue, for in hydrodynamical models of superwinds 
it is the merged SN ejecta that contains the ma jority of the 
energy of the wind (Stricklan d fc StevensH2000ri . and is the 
material most likely to escape into the IGM. If a significant 
fraction of the soft X-ray emission in superwinds comes from 
this material then future, higher-spectral resolution X-ray 
observatories, such as Astro-E2 and Constellation-X (which 



Conductive clouds in supernova- driven super-winds 21 



can measure the velocity of the hot gas) will be able to di- 
rectly constrain metal and energy ejection from starbursts 
into the IGM. If the hot gas is predominantly ambient ISM 
then these observatories will only probe a smaller fraction 
of the metal flow associated with superwinds. 

Observationally, the simplest diagnostic of the origin 
of the X-ray emitting gas would be to measure the absolute 
metal abundance. This has proved to be difficult, given a va- 
riety of systematic biases in the low to moderate resolution 
spectroscopy available with ROSAT, ASCA, Chandra and 
XMM-Newton. These in clude the artificial bl e nding of spec- 
trally distinct regions jPahlem et al.l 1 19981 : IWeaver et alJ 
120001 : iDahlem et al.l 12000^ . calibration uncertainties and 
over-sim plistic spectral model fits t o intrinsically complex 
spectra iStrickland fc Stevenslll998ri . The unphysically-low 
X-ray derived abundances found in th e earliest Chandra 
observations of supe r winds (see e.g . IWang et al.l l200li 
IStrickland et alJl2002l ; IXia et al"ll2002f) are almost certainly 
primarily due to small inaccuracies (< 2%) in the ACIS en- 
ergy scale available at that time. Re-analysis of that data 
with more modern calibrations leads to significantly larger 
best fit abundances that are statistically inconsistent with 
the previous results, although it is still effectively impossible 
to distinguish between Solar or super- Solar element abun- 
dances with respect t o hydrogen (e.g. iMartin et al.ll2002l 
IStrickland et, alJl2004h . In the near future observations with 
the much higher spectral resolution provided by the Astro- 
E2 calorimeter will not suffer from these problems, and so 
should provide strong constraints on the absolute elemental 
abundances in the nearest bright superwinds. 

We discussed the origin of the soft X-ray emission in 
our models of a wind/cloud interaction in Section f3.4l Here 
we restate those claims and develop our analysis further, by 
also considering cases where Z„ > Z c . 

Soft X-ray emission in the cloud models without ther- 
mal conduction (the NC models) generally arises in wind 
material within the bow shock (see Fig. [SJ, with a model- 
dependent contribution from the cloud surfaces. In the mod- 
els with conduction the soft X-ray luminosities are typically 
larger (by a factor between 2 and 30) than in the non- 
conductive models (Table |3J, and this increased luminosity 
is due to the increase in the hot gas density in the region 
between the cloud surface and the bow shock due to evapo- 
rated cloud material (see Fig. 1101 , 

Note that the Ovi luminosities of both classes of model 
are very similar. In these models the gas emitting Ovi typ- 
ically has a density of n ~ 5 x 10~ 2 cm -3 , and is always 
material originally from the cloud in both conductive and 
non-conductive models. As such it will have the metallic- 
ity of the unenriched cloud, i.e. low abundances if the cloud 
is a halo or high velocity clouds over-run by the wind, or 
the metallicity of the WIM of the host galaxy if it is being 
advected out of the disk within the wind. 

The results presented in Section 13.41 were made under 
the assumption that the metal abundance of the wind and 
cloud material were both Solar. As discussed in Section [5] 
the metal abundance of the wind material can be consid- 
erably higher than Solar, depending on the degree of wind 
mass-loading. In a situation where Z v > Zq wind mate- 
rial will make up a larger fraction of the X-ray emission, as 
the soft X-ray emissivity scales roughly in proportion to the 
metal abundance for gas in the temperature range T ~ 10 6 



- 10 7 K. Note that there must be little or no mass loading of 
the wind for its metal abundance to be high. Consequently 
the wind velocity and temperature will be higher and wind 
density lower than in a heavily mass-loaded case. Thus our 
Tl models could be considered as representing heavily mass 
loaded winds, and the T5 models as non-mass loaded or 
weakly mass loaded winds. 

To illustrate the influence of the wind metal abundance 
on the emission from the wind/cloud interaction, we recom- 
puted the X-ray luminosity as function of time in all mod- 
els assuming that the X-ray emissisivity of all wind mate- 
rial was either 4 or 10 times higher than for a Solar abun- 
dance plasma. In Fig. 1121 we show the ratio of this lumi- 
nosity Lx(Z) to the X-ray luminosities calculated assuming 
Z„ — Zq (which were shown in Fig.[5J as a function of time. 
Changing the abundance of the wind material will have neg- 
ligible effect on the Ovi emission, which is always dominated 
by cloud material. 

The large increase in X-ray luminosity seen in the 
T5LPNC models for high Z„ illustrates the dominance of 
wind material in this model (long dashed lines in Fig. 112b 
and b). Averaged over the entire length of the simulations 
wind fluid material accounts for 35% of the soft X-ray emis- 
sion when Z w = IZq. This fraction increases to 88% (95%) 
if Z„ — 4Zq (Zy, = IOZq). In contrast, with conduction 
active the wind material only accounts for on average 25% 
of the X-ray emission in model T5LP when Z„ — IZq. Al- 
though the increase in X-ray luminosity for this model at 
high Z w (see Fig. [TSJi and d) is significantly less than in 
model T5LPNC, wind material does dominate the emission 
(58% of the emission for Z„ = 4Zq , 77% of the emission for 
Z w = IOZq). 

All the Tl models are on average dominated by wind 
material, for all values of Z w , so the X-ray luminosity does 
increase strongly with increased Z w . 

The T5HP models are another interesting case. De- 
spite the very different dynamical state of the T5HPNC 
and T5HP models, the average fraction of the X-ray emis- 
sion from wind material in both models is relatively sim- 
ilar (although always slightly larger in the non-conductive 
NC models, representing the lack of any contribution due to 
strong evaporative loading of the bow shock region that is 
seen in the conductive model). For Z w = IZq wind material 
only produces 12% of the emission in the T5HP model, and 
14% of the emission in the T5HPNC. As argued previously 
a wind with conditions similar to those in the T5HP(NC) 
models can not be heavily mass-loaded, and hence will have 
a high metal abundance. Cloud material dominates in both 
of these models for Z w = 4Zq, but for Z w = IOZq wind 
material provides 57% of the emission in the T5HP model 
and 66% of the emission in model T5HPNC. 

Our discussion using broad energy band weighted lu- 
minosities and 1Z ratios is intended primarily as an aid to 
understanding, but some of the subtleties and second order 
effects are also worth considering. Not only does the rela- 
tive degree of enrichment of the wind (or Zsn/Zism) vary 
from element to element, but the degree of contribution of 
wind material (7?.) to the X-ray emission is also element 
dependent, given the different temperature of the conduc- 
tive interface compared to the bow shock. The emissivity 
of different elements and their ions are strongly tempera- 
ture dependent. In these simulations the conductive inter- 
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Figure 12. X-ray luminosities for high Z w normalized to Lx obtained with Z w = IZq. Black solid lines: TILP(NC); blue short dashed: 
TIHP(NC); red long dashed: T5LP(NC); green dot dashed: T5HP(NC). In all cases Z c = IZq. (a) Non-conductive models, assuming 
Z w = 4Zq. (b) Non-conductive models, assuming Z w = IOZq. (c) Conductive models, assuming Z w = 4Zq. (d) Conductive models, 
assuming Z w = IOZq. 



face has a lower temperature than the wind bow shock, 
and thus emission from lower-temperature ions (e.g. Ovn 
or Oviii) will be more heavily influenced by cloud material 
than higher temperature ions (e.g. Mgxi, Mgxn or Sixill). 
In some cases, particularly when conductive-interface and 
wind contributions are on average comparable, this effect 
would lead to peculiar X-ray-derived abundance patterns, 
with lower ionization state ions appearing to have lower 
abundances than high-ionization ions. It is conceivable that 
the peculiarly-low X-ray-derived ratio of Oxygen to o t her a - 
element abundances found in M82 bv lOriglia et al] (I2004T) 
might be caused by such an effect, although there are a few 
other observational biases that could also cause these odd 
results. 

The point that should be taken away is the gas probed 
by soft X-ray observations of superwinds is not necessary 
always un-enriched ambient gas even when thermal conduc- 
tion is active at the high Spitzer-conductivity levels. For 



some sets of conditions (and counter to naive expectation), 
in particular in powerful high velocity and high pressure 
winds, the ambient gas being conductively evaporated off 
the cloud does dominate the observable soft X-ray emission. 
In less powerful winds conductive evaporation is reduced, 
and in these case and in the case of no conduction, the soft 
X-ray emission arises in wind material within the bow shock 
around each cloud. The metal abundance of the wind ma- 
terial depends strongly on the degree of mass-loading of the 
volume-filling wind-fluid, which is a separate topic deserv- 
ing some additional investigation. However, this implies that 
if we can measure the metal abundance of the soft X-ray 
emitting gas it can be representative of the more tenuous 
volume-filling wind, and thus can be used to constrain the 
degree of large-scale wind mass loading. 
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4.4 Soft X-ray and Ovi emission and absorption 

4-4-1 Ovi absorption 

The mean Ovi column densities through a simulated sin- 
gle wind/cloud interaction (Table are typically 1.5 dex 
lower than the value of Now seen in starbursts observed 
with FUSE, which ra nge from 14.3 < logJVovi < 15.3 (see 
iHeckman et al.ll2002h . Increasing the cloud radius, or the 
metal abundance of the wind, is not likely to increase iVovi, 
for the reasons discussed in Section f3. 61 and Section f4. 31 

There is relatively little difference in mean iVovi be- 
tween the conductive and non-conductive models. While it 
is hard to assess to what degree numerical diffusion is arti- 
ficially increasing the amount of Ovi absorbing material in 
the non-conductive models, it is certainly true that Novi is 
not currently a strong discriminant between conductive and 
non-conductive models for wind/cloud interactions. 

The hotter and/or faster wind models produce larger 
mean iVovi per cloud. A high mean iVovi per cloud is advan- 
tageous, in that it reduces the number of clouds needed to 
match observations. Unfortunately, cloud stripping or evap- 
oration is more efficient in such models, so either more clouds 
would be required initially, or a continuous source of clouds 
would be required (e.g. via Kelvin-Helmholtz instabilities 
along the walls of the cavity). 

There would need to be ~ 30 clouds per line of sight if 
clouds alone are responsible for the Ovi absorption in super- 
winds. In practice the walls of the outflow will be another 
significant source of Ovi absorption. That there should be 
multiple clouds along any line of sight through a superwind 
is by no means surprising, given that the velocity width of 
observed optical and UV absorption lines is much greater 
than could be expected from any single structure. Neverthe- 
less, a model requiring very large numbers of clouds along 
any line of sight is neither parsimonious nor elegant. In sum- 
mary, the simulated Ovi column densities alone do not pro- 
vide strong constraints on allowable wind models when com- 
pared to the exi sting observation s. 

As shown in iHeckman et all J2002T) superwinds follow a 
distinct Ovi column density vs line width (b) relationship. 
We would expect that the column density and line width 
of a simulated Ovi profile from a single cloud should fall on 
the trends shown in Fig. 1 of that paper if a superposition of 
clouds is to create the net profile observed. Of the simulated 
Ovi profiles we investigated only h alf of them fell on the 
low Nov i, low b extrapolation of the IHeckman et all (I2002T) 
trend. The other simulated profiles had low Novi but high 
b values of ~ 200 km/s. This is puzzling, but it is too early 
to assess if this is a significant discrepancy. Our simulations 
only cover the first 1 Myr of a cloud/wind interaction, while 
real superwinds should contains clouds with a broad range 
of ages. 



4-4-8 Ovi and soft X-ray emission 

In Table |H] we summarize the published values of Ovi and 
soft X-ray fluxes in three starburst-driven superwinds (M82 
and NGC 3079 host strong superwinds, while NGC 4631 is a 
weaker starburst) and the marginal Ovi detection in the halo 
of a normal spiral galaxy believed to ha ve a galactic fountain 
(NGC 891), adapted from the work of iHoopes et alJ (I2003T) 



and lOtte et alJ i2003T) . We choose to concentrate our analy- 
sis on the ratio of Ovi to soft X-ray emission, and the ratio of 
the emitted FUV and X-ray power to the wind mechanical 
power, given that we can only simulate a single cloud. The 
30"-wide aperture of the FUSE observations samples a small 
fraction of the total region covered by soft X-ray emission 
in these objects, but this is considerably larger than volume 
we can simulate with high resolution. 

Note even with the limited observational data currently 
available it is clear that there is not a single universal 
Loyi/Lx ratio in superwinds: the upper limits on the ra- 
tio in M82 and N3079 are significantly lower than the values 
based on a statistically strong detection of Ovi in the NGC 
4631 wind. 

The value of the wind mechanical luminosity inter- 
cepted by the cloud, used in Table [21 is calculated over a 
radius of 10 pc, which is roughly the mean radius the com- 
pressed or crushed clouds assume. Note that the X-ray and 
Ovi luminosities are calculated using a larger radius of 50 
pc (in order to encompass the majority of the bow shock 
around the cloud). This mismatch in radii is reasonable, as 
the X-ray emission region is always larger than the cloud 
itself. If we were to use a 50 pc radius for both L mcc and the 
emitted radiation then the mean Lovi/imec and Lx/imec 
ratios would be lower by a factor 25. 

By using a radius equivalent to the effective cloud ra- 
dius in calculating L moc , and hence the Lovi/imoc, ix/imec 
and Ltot/imcc ratios, we are effectively assuming the cloud 
area covering factor (as seen by the wind) is of order unity. 
If the cloud covering factor is less than unity, then the mean 
radiative to mechanical energy ratios should be reduced pro- 
portionally. However, it should be noted that FUV spectra 
of starbursts taken with FUSE indicate that the area cover- 
ing fraction of neutral atomic cloud s (e.g. seen in Cll A103 6 
absorption) is indeed of order unity iHeckman et alJboOlal) . 
in contrast to the lower area cove ring fraction of (pre sum- 
ably denser) molecular gas clouds llHoopes et al.ll2004l) . 

Table [3] shows the mean fraction of wind mechanical 
energy converted in X-ray or FUV radiation by a single cloud 
of initial radius 15 pc, along with the mean Ovi to soft X- 
ray flux ratio. As Ovi emission accounts for ~ 30% of the 
total emissivity of gas at T ~ 3 x 10 5 K, we can combine the 
Ovi luminosities with the soft X-ray luminosities to gain an 
estimate the of total radiative energy loss from the coronal 
and hot phases. 

We find that a low fraction of the wind mechanical en- 
ergy intercepted by a single cloud is lost to radiation. The 
total fraction of wind energy radiated per cloud found in 
the models assuming Z w = IZq lies in the range ~ 0.25 
- 1%, and is relatively similar between the conductive and 
non-conductive models. For higher values of Z w the ratio 
ix/^mec per cloud is still < 1%, except in the T1LP and 
T1LPNC models where it reaches 2 - 3%. As mentioned 
previously a somehwat larger fraction of wind energy ends 
up as cloud kinetic energy. 

The simulations show that a superwind can flow around 
a cloud and thus affect multiple clouds along any particular 
path without any substantial radiative loss of wind energy, 
except in the case of the TILP(NC) models where radia- 
tive losses may violate observational limits. Summed over 
all clouds the Ltot/£mec ratio should not significantly ex- 
ceed estimates of the X-ray and FUV luminosity to wind 
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mechanical power in real superwinds. Usi ng FUSE obser- 
vations of the dwarf starburst NGC 1705, Hcc kman et alJ 
(2001b) limit to the coronal cooling rate to be ~ 5/e% of 
the mechanical luminosity based on the observed Ovi ab- 
sorption line strength. Soft X-ray observations of superwinds 
give a ratio of ~ 1/e - 3/e% fo r the X-ray to mechani- 
cal power ratio dStricklandl 12004 1 , The major uncertainty 
in these observationally-derived ratios is the unknown effi- 
ciency of supernovae energy thermalization e. 

The lower pressure wind models (LP models) produce 
more emission per unit mechanical power than the high pres- 
sure models, although the magnitude of Lovi or Lx is larger 
per cloud in the high pressure models. The lower tempera- 
ture wind models (all Tl models) also produce more emis- 
sion per unit Z/ mec than the higher temperature wind models, 
but this is mainly because of the high wind density in these 
models. 

In accord with the observational evidence, coronal phase 
cooling (i.e. emission predominantly in the FUV) dominates 
over X-ray cooling from hotter gas by a factor of at least 
a few. The strongest difference between the models consists 
in the X-ray cooling efficiencies. With the one exception of 
model T1LPNC, the Z/x/L mec ratios are an order of magni- 
tude lower in non-conductive models compared to the mod- 
els with conduction. 

Another significant difference between the conductive 
and non-conductive clouds models is in the Ovi to X-ray 
luminosity ratio. In our opinion this ratio places the most 
interesting constraints on whether wind/cloud interaction 
models can reproduce observed superwind properties. 

For Z„ — Zq the majority of the non-conductive models 
have Loyi/Lx ratios that are orders of magnitude larger 
than the largest values observed to date (compare Table[3]to 
Table El . The low pressure non-conductive model T1LPNC 
has Lovi/Lx = 1.07, but this value is still much larger than 
the upper limits on the Lovi/ix ratio for M82 and NGC 
3079. 

The conductive models span a much smaller range in 
Lovi/Lx, from 0.2 to 2.7. Again, the low temperature low 
ram pressure models have the lowest Z/ovi/Lx ratios, closer 
to the observed range, if still only matching the values ob- 
served in NGC 4631. 

The fundamental problem with the non-conductive 
models are the low X-ray luminosities. The Ovi luminosi- 
ties are comparable in both conductive and non-conductive 
models, so it is the higher soft X-ray luminosities generated 
by cloud material in the conductive interface that leads to 
roughly acceptable Lovi/Lx ratios. It may be the case that 
the typical Lovi/Lx ratio in superwinds is < 0.1, in which 
all of our wind/cloud interaction models have trouble pro- 
ducing a realistic Lovi/ix ratio if Z v ~ Z c . 

Larger clouds are more effective at intercepting wind en- 
ergy and converting it into radiation. The mechanical power 
intercepted scales as R%, but the X-ray emission scales as 
roughly R% over the small range of radii we have investi- 
gated (see Section [3. 611 . The Ovi luminosity does not grow 
with increased radius as rapidly as the X-ray luminosity, 
so the Loyi/Lx ratio per cloud will drop as for clouds of 
larger radius. This effect is probably not significant enough 
to reconcile the low Z v non-conductive models with ob- 
served Lovi/Lx ratios, but should be considered for the 
conductive models. In the case of models with thermal con- 



duction, matching the Lqyi/Lx ratio in M82 and NGC 3079 
is possible if the average cloud has an initial radius only a 
few times larger than 15 pc. Larger cloud radii also have 
the added advantage of longer cloud lifetimes, as previously 
discussed. 

Even values of Z„ = IOZq are not sufficient to pro- 
duce acceptably low Lovi/ix for the non-conductive mod- 
els (see Table [3J|. However, for the conductive models a 
high Z w value does produce Lovi/Lx values similar to the 
Lovi/Lx < 0.1 upper limit in M82 and NGC 3079. Thus 
even in the case of Z v ^> Z c the Ovi to X-ray ratio appears 
to provide a strong discriminant between non-conductive 
and conductive cloud models for superwinds, with conduc- 
tive models producing values in the range of those observed. 

Furthermore, it appears that we need not require con- 
duction at full Spitzer levels to attain acceptably low 
Lovi/Lx ratios. As shown in Table a reduction in the 
value of the coefficient of thermal conductivity by a factor 
/ = 5 barely altered the time-averaged Lovi/Lx ratio. Even 
for a reduction of / = 25 the increase in Lovi/£x is sur- 
prisingly small. Conduction is still dynamically important in 
these models in stabilizing the cloud against fragmentation, 
so it is possible that with weak to moderate conduction we 
can produce longer- lived clouds with suitably low Lovi/ix 
ratios. Further simulations tracking clouds over significantly 
longer timescale will be necessary to test this hypothesis. 

For the set of parameters investigated in this paper the 
models T1HP and T5LP have the most acceptable X-ray 
and Ovi emission properties, in that the Ltot/imec ratio 
per cloud is not too large (even in cases where Z w ~ IOZq) 
and the Lqyi/Lx ratio is also low. 



5 SUMMARY 

We performed a series of 2-dimensional hydrodynamical sim- 
ulations with a resolution 0.1 pc of a super-or-Hyper-sonic 
starburst-driven superwind interacting with an initially sta- 
tionary cool dense cloud of radius R c = 15 pc. Such interac- 
tions are one of a class of wind/cloud interactions known to 
occur within superwinds. Wind/cloud interactions may play 
a significant role in shaping their observable properties, in 
particular in the FUV and X-ray energy bands, but have 
not been explored numerically in any detail. More general 
analytical models of ten applied to superwinds, such as the 
I Weaver et al] il977t) wind-blown bubble model, have trou- 
ble explaining observed Ovi absorption line column densities 
and Ovi to X-ray emission flux ratios. 

We considered a range of models, with low and high 
wind ram pressure (both pressures are much higher than 
conventional ISM pressures), low and high wind to cloud 
density ratios, and the presence or absence of thermal con- 
duction. The resolution of the simulations is high enough 
that the region of the conductive interface in which the 
strong coolant Ovi is produced is numerically resolved. A 
numerical tracer field was used to distinguish between gas 
originally part of the cloud and gas from the wind. We fol- 
lowed the evolution of the wind/cloud interaction for times 
up to 1.5 Myr, at which point the cloud has been carried off 
the numerical grid by the wind. In addition to considering 
the dynamics of the wind/cloud interactions, we have also 
modeled their soft X-ray emission and FUV Ovi emission 
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and absorption properties and compared them to existing 
observational data. 

There are significant dynamical differences between 
wind/cloud interactions in the presence of thermal conduc- 
tion compared to those without conduction. In the absence 
of conduction the clouds are fragmented in a few dynam- 
ical (cloud-crushing) time-scales t cc = X^^c/iAv (where 
X — Pc/pw), as ex pected from earl ier of clouds over-run 
by SN blast waves llKlein et alJll994|) . Relatively little cool 
cloud mass is converted into hot gas by shocks or mixing, 
but the cloud rapidly ceases to exist as a recognizable co- 
herent object due to the hydrodynamical instability-driven 
fragmentation. 

In contrast, thermal conduction acts to stabilize the 
cloud by inhibiting the growth of Kelvin-Helmholtz and 
Rayleigh- Taylor instabilities (as had previously been found 
in the case sub-sonic motion of M olecular clouds in a hot 
medium bv lVieser fc Henslerll200(iF) . Furthermore, the con- 
ductive heat flux into the cloud's surface generates a shock 
wave that compresses the cloud and modifies its internal 
density structure significantly and reduces its conductive 
m ass loss rate (as previously di scovered for stationary clouds 
bv lFerrara fc Shchekinovll993l) . These dynamical effects oc- 
cur even in models in which the coefficient of thermal con- 
ductivity has been reduced by factors of up to 25 from the 
Spitzer value. If thermal conduction is active, then clouds 
embedded within a superwind will exist as a single coherent 
structure for longer time scales than they otherwise would, 
but at the cost of losing mass to evaporative heating. The 
strength of these conductive effects is proportional to the 
dimensionless parameter ao (Equation|5J. This parameter is 
high for low density winds, small clouds, and high wind tem- 
perature, especially in the case of supersonic winds where 
the temperature of material in the bow shock around the 
cloud is greater than the wind temperature. 

In general, the higher the wind ram pressure or velocity 
is, the higher the mass loss rate of a cloud will be, due to 
either hydrodynamical fragmentation or conductive evapo- 
ration. 

With or without the presence of thermal conduction, 
Ovi emitting and absorbing material primarily arises in 
cloud material at the periphery of the cloud and any cloud 
fragments. In contrast, the origin of the soft X-ray emitting 
material depends on the strength of conductive effects and 
the relative metal abundances of the wind and cloud ma- 
terial. If conductive effects are weak or absent, and/or the 
metal abundance of the wind Z w is highly enriched com- 
pared to the cloud Z c , then the observed soft X-ray emis- 
sion comes from wind material in the bow shock around 
the cloud. Thus even if the undisturbed superwind is rela- 
tively tenuous and thus a weak X-ray emitter, when com- 
pressed around an obstacle such as a cloud its emission is 
enhanced and hence the wind abundance pattern could be 
probed by soft X-ray observations. For strong conduction, 
or in the case of weaker conduction with Z w ~ Z c (e.g. a 
heavily mass-loaded superwind) the X-ray emission again 
arises within the bow shock region, but conductively evap- 
orated cloud material mixed into this region dominates the 
emissivity. In the highest wind ram pressure and wind tem- 
perature models with conduction or without conduction it is 
still cloud material that dominates the X-ray emission, even 



if the metal abundance of the wind is 10 times greater than 
that of the cloud. 

Thus whether soft X-ray observations probe wind or 
cloud material in a wind/cloud interaction depends on sev- 
eral factors: primarily the strength of conduction (conduc- 
tion favouring cloud material) and the relative abundances 
of wind with respect to the cloud (if the wind is strongly en- 
riched then its influence increases) . Different material can in 
principle been seen at different X-ray energies, with low ion- 
ization species such as Oxygen being more strongly weighted 
toward probing the cooler conductive interface than high 
ionization species such as Mg or Si (which would tend to 
probe the hotter regions of the bow shock) . This complexity 
may contribute to some of the uncertainties and peculiar- 
ities found in c urrent X-ray-deriv ed metal abundances in 
starbursts (e.g. lOriglia et al.ll2004) . 

The Ovi absorption line column density produced per 
cloud is relatively similar between conductive and non- 
conductive models, ranging from 12.9 < log iVovi < 13.5. 
In all cases the mean iVovi value found is substantially 
less than the observed total Ovi column densities in star- 
bursts with winds, so that in the wind/cloud interaction 
model requires many clouds per line of sight (typically 10 - 
30 clouds). Models with stronger hydrodynamical stripping 
or conductive evaporation produce higher mean Novr, but 
their shorter cloud lifetimes are disadvantageous for accumu- 
lating larger numbers of clouds, and maintaining cool cloud 
gas populations in winds for the ~ 10 Myr time periods that 
appear to be implied by observations. 

In terms of emission we considered two primary obser- 
vational constraints. The strongest observation constrain is 
the ratio between Ovi and soft X-ray emission Lovi/^x, 
whic h has the va lue ~ 1 in NGC 4631 (this work, see also 
lOtte et aljl2003l) . and is definitely much lo wer (< 0.1) in 
the s tronger M82 and NGC 3079 superwinds llHoopes et alJ 
I2003L Hoopes et al, in preparation). The second is that 
observations place limits on the fraction of wind mechan- 
ical power that is radiated in the FUV and X-ray bands, 
which are relatively low (e.g. < 5/e%) but made somewhat 
uncertain by our lack of robust knowledge of the starburst 
SN thermalization efficiency e. We found that except for the 
low ram pressure low wind temperature TILP(NC) models, 
< 1% of the wind mechanical energy is lost to radiation per 
cloud, even for wind metal abundances as high as WZq. 
Typically between 2 - 6% of the wind mechanical energy 
ends up as kinetic energy of the cold cloud within the first 
0.5 Myr. 

If we consider both possible cases, of both Solar and 
super-Solar abundance winds, we can draw the following 
conclusions. Qualitatively, the high wind-ram pressure mod- 
els or the low ram pressure high wind temperature model are 
somewhat better than the low pressure, low wind tempera- 
ture, models at matching the available constraints on X-ray 
emission and Ovi emission and absorption. If Z v ~ Z c then 
the non-conductive models are ruled out by their very high 
Lavi/Lx ratios (due to the relative lack of X-ray emission), 
while conductive models can match the value of ~ 1 ob- 
served in NGC 4631. Even if Z„ = 10Z e then the Lovi/L x 
ratio in most of the non-conductive models remains higher 
than the value observed in NGC 4631. For high Z w the con- 
ductive models produce Lovi/Lx as low as ~ 0.1, similar to 
the upper limits on the ratio in M82 and NGC 3079. 
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Thus some, but not all, of our current wind/cloud in- 
teraction models can match the range of Lovi/ix values 
observed in the starbursts NGC 4631, M82 and NGC 3079, 
if conduction is active, and/or with high wind metallicity 
(i.e. little wind mass loading). The low wind ram pressure, 
low wind temperature models (T1LP/NC) radiate too large 
a fraction of the wind's mechanical energy per cloud to allow 
large numbers of clouds along any line of sight, which may 
rule them out as a plausible model. A smaller fraction of the 
superwind's energy is radiated per cloud in the high pres- 
sure models, but the total luminosity per cloud is slightly 
higher. This is advantageous, as fewer clouds are required 
instantaneously, but the ultimate issue of cloud lifetime and 
replacement still remains. Larger clouds or lower thermal 
conductivity are plausible ways in which cloud lifetimes can 
be extended to the ~ 10 Myr time scales required by obser- 
vations. 

Wind/cloud interactions in superwinds have been dis- 
cussed qualitatively for a long t ime, stretching at l east a s 
far back as the seminal paper of IChevalier fc Cleed Il985ft . 
The conditions that appear to favor long cloud life times 
are not necessarily those that appear to produce the best 
match to FUV and X-ray emission and absorption proper- 
ties, but it does appear that these competing requirements 
can be balanced. More work is required to follow clouds for 
longer periods, and to better investigate cloud acceleration 
and thermal conduction at sub-Spitzer levels. 

That there is some level of diversity in the physical 
conditions within different starburst-driven superwinds and 
thus in their emitted radiation is not totally unexpected. In 
principle this can be advantageous if we are to disentangle 
the multiple fluids variables that shape a superwind and un- 
derstand their physics. Observing more superwinds in detail, 
especially with FUSE, will be scientifically profitable. 

Nevertheless, we find it very encouraging that some of 
these wind/cloud models can quantitatively match some of 
the observed properties of winds. Furthermore, these simu- 
lations have been very enlightening with regard to the both 
the role conduction may play in winds (a topic that has seen 
very little direct work despite its obvious potential impor- 
tance), and to what controls which material (cloud or wind) 
is probed by X-ray observations. This latter issue is apposite 
given the immanent launch of the very high spectral resolu- 
tion X-ray micro-calorimeter on Astro-E2. With direct ob- 
servational probes of the elemental composition and possibly 
the kinematics of the soft X-ray-emitting gas in superwinds 
it will be possible to constrain the physical processes that 
shape emission from these outflows. 
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